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Abstract 

I review the treatment of high-energy QCD in Minkowski space, with an emphasis on 
factorization theorems as extensions of the operator product expansion. I discuss how 
the factorization properties of high-energy cross sections and amplitudes lead to evolution 
equations that resum large logarithms for two-scale problems. 



1 Based on seven lectures at the Theoretical Advanced Study Institute, QCD and Beyond, Boulder, Colorado, 
June 1995. 
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1 Introduction: Fundamental Issues and Problems 

Many of the fundamental properties of field theory, such as the operator product expansion, are 
best developed in Euclidean space. The very nature of space-time, however, dictates that an 
understanding of field theory directly in Minkowski space is also indispensable. I hope these lec- 
tures will arm the reader with general insights into how perturbative quantum chromodynamics 
(QCD) manifests itself at high energy in Minkowski space. I will discuss the calculational chal- 
lenges that result from singularities on the light cone, and develop some of the methods that 
we currently possess to deal with them. 

The unifying thread that runs through these lectures is factorization [|TJ, the systematic 
separation of dynamics associated with short and long distance scales [0, |3| . This is a recurring 
theme in modern theoretical physics, which can also be found as a central idea, sometimes under 
the names effective actions or field theories, in other lecture series in this school. One aim of 
these lectures is to identify quantities in QCD that are genuinely short-distance dominated, 
which will lead us to the concept of infrared safety, and the ubiquity of jet cross sections @]. 
Another is to show how one of the great phenomenological successes of high energy physics, 
the parton model || [|, emerges as a consequence of the factorization properties of QCD in 
Minkowski space, which lead as well to its systematic improvement, including evolution [[/J. 

Resummation refers below to the summation of enhancements (usually logarithmic) in 
ratios of kinematic variables, such as energy and momentum transfer, to all orders in field- 
theoretic perturbation theory. Rather than try to review all recent progress in this large and 
growing field, I will emphasize two representative and classic examples, Sudakov [0] resumma- 
tion in e + e~ annihilation and BFKL || resummation in deeply inelastic scattering, stressing 
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how they may be regarded as consequences of the underlying factorization properties of field 
theory. 

These lectures are concerned primarily with the theoretical foundation of perturbative 
QCD. Of necessity, much introductory material involves quite general properties of quantum 
field theory. Much of this material may be found in field theory textbooks |TD|, [TTJ , as well 
as in the books by Yndurain and Muta on QCD [Q. although some developments are given 
below which may be less widely familiar (Tl|, [13] . In particular, I have tried to outline a general 
analysis of long-distance behavior in perturbation theory, based upon the analytic structure 
of Feynman diagrams |14| and on an infrared power-counting procedure. The reader will find 
many points of contact, but many differences in emphasis, with my TASI 91 lectures ||15|| . 

Many of the important and innovative calculational techniques and the phenomenological 
analyses that realize the program described here, as well as other equally important aspects 
of the theory, are treated in other lectures at this school. Very useful recent reviews that 
treat these topics in a more directly phenomenological manner include the TASI 94 lectures of 
Ellis |T(|, and a "Handbook" by the CTEQ Collaboration [jT7|| . As a review of perturbative 
QCD, the collection of monographs edited by Mueller [plj is extermely helpful, and introduces 



and reviews subjects, particularly elastic scattering [19| and QCD coherence |20|, [21]], that are 
closely connected to the discussion in these lectures. Finally, for a theoretical introduction from 
a complementary point of view, see the book of Dokshitser, Khoze, Mueller and Troian p2]| . 

So much of the terminology of QCD is intertwined with its sources in experiment and 
theory, that it seems appropriate to begin with a very brief review of the strong interaction 
physics that led to QCD. This will be followed by an introduction to deeply inelastic scattering, 
the pivotal experiments whose outcome made it possible to identify quantum chromodynamics 
as a promising theory of the strong interactions over twenty years ago. The balance of these 
lectures will discuss how this promise has been realized. 



1.1 A Prehistory of QCD 

It is plausible to identify the birth of strong interaction physics with the discovery of the 
neutron, which along with its beta decay, signaled the existence of two new interactions, the 
strong and weak. From this slender thread, physicists suspended the four-fermion theory of 
weak interactions, and the concepts of (strong) isospin, the S-matrix, and the Yukawa theory, 
in which the strong interaction is mediated by the exchange of scalar particles, now identified 
as pions, between nucleons. Beyond this, through the thirties and forties into the fifties, strong 
interaction physics remained for the most part nuclear physics, with many well-known successes 
and consequences pjfl . 

With the discovery of the pion and then other mesons, around 1950, and the development 
of high energy accelerators, the substructure of nucleons became for the first time an object 
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of study. The earliest information was on the excited states of mesons and baryons, and 
the discovery of the A in pion-nucleon scattering ushered in a bewildering array of resonances. 
Order was imposed on this chaos through symmetry principles, leading eventually to the concept 
of quarks as the building blocks of hadrons; three quarks for baryons and a quark-antiquark pair 
for mesons. This "constituent" quark model still successfully describes most of the qualitative 
features of baryon spectroscopy [£4], 25]. It was thus natural to try and "see" the quarks in 
experiments with momentum transfers large enough to resolve the internal structure of the 
nucleon, and to explore the possibility that forces between quarks are mediated by a field, 
dubbed, perhaps for want of a better name, the gluon. 

At this time it was by no means obvious, or universally recognized, that such a picture of 
strong interactions would or could succeed. Indeed, it was a widely held view that the idea of 
elementary particles was inappropriate for the strong interactions altogether, and that theories, 
particularly those involving perturbative methods, could serve only as guides to suggest the 
mathematical properties of the S-matrix [[U|] . 

The deeply inelastic scattering experiments, which we shall discuss shortly, changed all 
that. They showed unequivocally that the proton possessed charged substructure of a spatial 
size much smaller than the proton itself. Indeed, the experiments also suggested spin one-half 
for these particles. At the same time, it was recognized that the constituent-quark model 
seemed to require a new quantum number for the quarks, "color" . Color, although originally 
introduced to solve the problem of Fermi statistics for the spin-1/2 quarks, provided a natural 
set of currents to which the gluons might couple. A three-color model of quarks has a global 
SU(3) symmetry, with currents reflecting its group structure. This suggested a local nonabelian 



gauge theory of the type originally introduced by Yang and Mills many years before [26 



In a relatively few months this theory was recognized to possess a number of important 
properties. At lowest order, it produces attractive forces in the three-quark and quark-antiquark 
systems, and repulsive forces for a quark-quark system. It automatically incorporated the well- 
established successes of current algebra |. Perhaps its characteristic and most crucial property, 
however, is asymptotic freedom, according to which its coupling decreases with decreases in the 
distance scale over which it is measured p8| . The need for asymptotic freedom was signalled 
by the very experiments designed to detect point-like structure in the nucleon. These are the 
deeply inelastic scattering (DIS) experiments, initiated at SLAC in the late sixties, to which 
we now turn. 



2 In essence, current algebra follows from the assumption that the electromagnetic and weak interactions 
couple to hadrons via point-like operators (the currents) which are then assumed to obey commutation relations 
consistent with the symmetries of these interactions (the algebra). In QCD, this algebra follows automatically 
by identifying the currents with operators qi'y^qj, with qi quarks of differing flavors. 



5 



1.2 Deeply inelastic scattering and the parton model 



In deeply inelastic (deep-inelastic) scattering, a massive hadronic state X of invariant mass 
M\ ^> m 2 N is produced by the scattering of a lepton (for instance, an electron) on a nucleon 
(or other hadron), 

e{k) + N(p) -> e(fc') + X hadronic . (1) 
This process in illustrated in fig. 1. 




X 



Figure 1: Deeply inelastic scattering. 

Kinematics. Because the lepton interacts with the nucleon only through the exchange 
of a photon, W or Z, with relatively small electroweak corrections, the cross section for this 
process factors into leptonic and hadronic tensors, 

d z h' 1 

rfff = aW r(M) ^ (M) ' (2) 
where we have taken the example of photon exchange. The leptonic tensor L^ u is known from 
the electroweak Lagrangian, while the hadronic tensor may be expressed in terms of matrix 
elements of the electroweak currents to which the vector bosons couple, 

IT ee ^-tr[^'Y} 

W, u = -L £ Y,< n M\M q )\ x >< x \Uo)\n(p,<t)> 



8tt . 

spins a X 



X{2^5\p x -q-p). (3) 



The expression for L^ v is elementary, and the expression for is quite general, depending 
only on the assumption that the electroweak interactions couple to the hadron via local currents 

B 

3 This assumption, we have observed above, had considerable experimental and theoretical support by the 
late sixties. Nevertheless, it is itself a nontrivial assertion. 
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The matrix elements in W^ v include the strong interactions only, and they hence satisfy 
symmetry properties of the strong interactions. For instance, both electromagnetic and strong 
interactions enjoy invariance under parity and time reversal. This leads to a symmetric 
when the hadron is unpolarized, while time reversal invariance of the strong interactions leads 
to a real hadronic tensor, 

KT ] = W »7 ] ( § P in " averaged), W, v = W% . (4) 

Along with Lorentz invariance, and electromagnetic current conservation, 

q»W, u = 0, (5) 

these constraints may easily be used to show that the sixteen components of the hadronic tensor 
are determined by two independent structure functions, W\ and W 2 , 



The Ws are functions of Q 2 and the dimensionless ratio, 

x = = . (7) 

2p ■ q 2p ■ q 

For reasons which will become clear in a moment, it is also convenient to introduce dimensionless 
structure functions, 

Fi = Wi, F 2 =p-qW 2 . (8) 

It is worth noting that in the literature many definitions are given for the W s (usually differing 
by factors of target mass), while definitions of the F's are much more standardized. 

Scaling. The striking result of the early deeply inelastic scattering experiments was that, 
for Q 2 > 1 GeV 2 , the structure functions F(x,Q 2 ) become functions of x only, nearly inde- 
pendent of Q 2 . This property is called scaling |2£], |50fl . 111 hs importance, it is (and should 
be) compared to the experiments of Rutherford and collaborators, who discovered the atomic 
nucleus in the wide-angle scattering of alpha particles. In deeply inelastic scattering, the wide- 
angle scattering of the electron serves to detect point-like structure in the nucleon though 
scaling. If electic charge were uniformly distributed within the nucleon, we would expect wide- 
angle scattering to be very rare, giving structure functions that decrease rapidly with Q 2 . To 
see how point charges produce scaling, we shall review the parton model. 

The parton model. In the parton model ]5], |(| , we imagine the proton, or any other hadron, 
to be made of point-like constituents, the partons, through which it couples to the electroweak 
interactions, and eventually to the strong interactions as well. Deeply inelastic scattering will 
give information on the spin of these partons, which, in anticipation, we take to be one-half. 
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The fundamental relation of the parton model for deeply inelastic scattering may be 
written 

d^ N \p, q) = £ f da Bo JW(fr, q^fMO , (9) 
/ 

where da^ N \p, q) is the inclusive cross section for nucleon-electron scattering, while daBom(£,P, 
q) is the lowest-order, elastic parton-electron cross section, with the parton's momentum given 
by £,P, £, between zero and one. The functions 0//jv(O are parton distributions, which describe 
the probability of finding a parton (of "flavor" /) in the "target" hadron N. There is assumed 
to be no interference, either between different flavors or between different fractions £. This 
"incoherence" is a hallmark of the parton model and its extension to QCD. An important 
consequence is that the parton distributions are universal, in the sense that, since they describe 
processes that do not interfere with the hard scattering, they are the same for all inclusive 
hard scattering processes, not only for electromagnetic DIS. The most obvious extension is to 
neutrino DIS, but there are many others. Eq. (||) is illustrated by fig. 2, in which the hadronic 
interaction of fig. 1 is broken up into a parton distribution and parton-electron scattering. 




Figure 2: Deeply inelastic scattering in the parton model. 

Scaling and spin in the parton model. Both the inclusive hadronic cross section and the 
elastic (exclusive) partonic cross section may be factored into leptonic and hadronic tensors as 
in eq. (0). The hadronic tensor at Born level is given explicitly by 

w,J f) = i / Tr^-QAy'^ (2tt) 4 5V - - ?) , (io) 

with Qf the electromagnetic charge of quark flavor /. Substituting fllCf ) into ([5]), and using the 
definitions (||) and @, we find the following simple relations between hadronic and partonic 
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structure functions, 

f J ° 

Fi {N \x) = $ F^ f \x/Oh/N(0, (11) 

from which we can read off the partonic structure functions, 

2F^\z) = F 2 ^(z)=Q f 2 5(l-z). (12) 

Note the minor difference in the powers of £ in the convolutions (|TTD , which is a reflection of 
the tensor structure of W^ v . 

Turning again to the general relation eq. (H), we derive nucleon structure functions in the 
parton model 

F2 {N \x) =Y,QMfMx) = 2xF 1 (N \x) . (13) 
/ 

The relation F 2 = 2xFi is a direct consequence of the spin of the partons. This result [j3l[l , 
known as the "Callan-Gross relation" is important evidence that the partons detected in deeply 
inelastic scattering are indeed the quarks of hadron spectroscopy. 

A crucial feature of eq. ([13]) is that the structure functions Fi are independent of Q 2 . 
This is the scaling of experiment, and is a direct result of the assumed point-like nature of the 
partons. In summary, the parton model serves both to explain scaling and to identify partons 
as quarks. 

Heuristic justification and requirements on a field theory. The ultimate justification of 
the parton model must be in quantum field theory, but many of its features can be understood 
on the basis of compelling heuristic argurments. Fig. |3| illustrates the time development of 
a deeply inelastic scattering event. Fig. |^a shows the system before the collision, as seen in 
the center of mass. The isolated electron approaches from the left, and from the right the 
nucleon, of momentum p. The nucleon is pictured as a set of partons, spread out more-or-less 
evenly in the transverse direction. The nucleon is highly Lorentz contracted in the longitudinal 
direction. We must first justify why the nucleon may be treated as a "pure" state with a definite 
number of partons. For this, we imagine that the wave function of these partons is formed by 
interactions which occur on time scales of the order of 1 GeV -1 . Such scales are long compared 
with the time it takes the electron to traverse the nucleon, and hence remain uncorrelated with 
hard processes that the electron may initiate in that time. During that short time, each of 
the partons may be thought of as possessing a definite fraction of the nucleon's momentum, 
denoted £p in the figure. 

To initiate a hard scattering, the electron should pass very close to one of the partons, i. e. 
at 0(Q~ l ), close enough so that they may exchange a photon with an invariant mass Q 2 . Such 
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Figure 3: Heuristic picture of parton model DIS. 



a collision is shown in fig. |3|b. Actually, in most cases the electron will miss the highly localized 
partons altogether, and if it does come close to one, it is highly unlikely to hit another. Indeed, 
the probability for encountering n partons at a distance of Q^ 1 behaves as 



partons 



(14) 



with Ro the nucleon radius. Multi-parton events are thus power-suppressed in the parton model, 
a feature that will turn out to have a natural field-theoretic analog. Fig. |3]c shows the results 
of the collision. Long-distance forces produce the observed hadrons from the fragments of the 
struck nucleon. These forces, however, are again on the scale of hadronic internal interactions, 
and hence cannot interfere with the short- distance scattering of fig. |||b. In a rough summary, 
deeply inelastic scattering shows incoherence because initial-state interactions (which bind the 
proton) are too early, and final-state interactions (which produce final-state hadrons) are too 
late relative to the short time scale of the hard scattering. 

The above arguments rely on very general features of relativity and quantum mechanics. 
To justify their application within a specific quantum field theory, however, is highly nontrivial. 
The particular challenge concerns the separation of short and long time scales. We have assumed 
that at short scales the "strong" interactions, which evidently produce confinement, are weak 
enough to be neglected compared to the single, electroweak interaction. As we have seen in 
Collins' lectures at this school, interaction strengths that change with distance scales are the 
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rule in quantum field theory. It is clear, however, that the successes of the parton model, 
and their consistency with long-distance strong interactions, requires an asymptotically free 
theory. Thus, we are ready to review elements of the asymptotically free theory of the strong 
interactions, quantum chromodynamics. 

1.3 QCD, asymptotic freedom and the running mass 

Quantum chromodynamics is specified by its Lagrange density, as a nonabelian SU(3) gauge 
theory for Dirac spinors. Including gauge fixing and ghost terms the density for an SU(N) 
theory may be written as 

f 

- \tr{F, v Fn-\{ri-A a f 

+ vM d ^ S ad - gC abd A b »)c d . (15) 

Here the (Dirac spinor) quark fields are qj of mass Mf (/ = 1, . . . n/ labels flavor) with color 
index i, where i — 1 . . . N for an SU(N) theory. The gluon fields are A%, a = 1 . . . N 2 — 1, and 
ghost fields c a (and antighost c a ), a = 1 . . . N 2 — 1. The are SU(N) generators in the 
iV-dimensional defining representation. In matrix notation, the nonabelaian "field strengths" 
Fp, = F^, a T^ F) are given by 

pv» = q» A v _ d v A » + a v \ , (16) 

where the commutator is in terms of the color matrices. A number of useful identities and 
other properties of the color generators are listed in Appendix A. 

In the gauge-fixing term (\/2)(rj ■ A) 2 , rj is typically chosen as the gradient d (covariant 
gauge) or as a fixed vector n. For the latter choice, it is sometimes convenient to distinguish 
the possibilities n 2 > (temporal gauge), n 2 < (axial) and n 2 = (light-cone). In the limit 
A — > oo, the condition r] ■ A = is enforced, and the gluon decouples from the ghost. 

The diagrammatic rules for perturbation theory that follow from Eq. ([To]) are derived by 
standard techniques |TD|, [TT], |T^] . Their most characteristic feature is the self- interaction of the 
massless gluons and, in covariant gauges, the necessity of ghost fields. 

The theory, of course, must be renormalized, and the strength of its coupling therefore 
depends upon ( "runs with" ) the momentum scale at which it is defined, according to 

q^9W \ 90 = P WJ) , (17) 

where the bare coupling go held fixed. The beta function /3(g) is determined to one loop from 
the coupling renormalization constant Z g by 
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(91n/i 



' - ln(/i 2 /M 2 ) J " 7 



5tt 2 



9 In// 
16vr 2 V 3 3 

where M 2 is an ultraviolet cutoff, and where in the final line we have used g = go to lowest 



— N — ~ rif I , (1* 



order. The solution to fllTD is then the familiar 



0>) = " 

1 + W°2ln^ 

IfW 2 

(19) 



6 2 ln(/i 2 /A 2 ) ' 

with 6 2 = 11 — 2rif/3 for QCD and where in the second form we define 

A = ^-s^fe^o) ) (20) 

independent of /i - Evidently, the coupling g(fi), or equivalently "alpha-strong", 

decreases for large momenta, or short distances (with corrections due to higher terms in (5(g)). 
At the same time, as the momentum scale decreases, and distances correspondingly increase, 
the perturbative coupling grows. This is just what we needed for the parton model. It is, 
however, still a long way from this observation to a phenomenology of deeply inelastic and 
other hard-scattering processes. Nevertheless, it is with asymptotic freedom that everything 
begins. Without it, there is no natural explanation in field theory of the successes of the parton 
model. 

The last item in this synopsis of QCD and asymptotic freedom is a brief discussion of quark 
masses. Unlike quantum electrodynamics, where the physical mass of the electron is directly 
observable, the masses of the quarks must remain to a large extent theoretical constructs, which 
may be determined from experiment once a scheme for doing so is defined, but for which there 
is no unique scheme. In this sense, quark masses are much like the gauge coupling itself, and 
for many purposes it is useful to define a running quark mass, 



m(fi 2 ) = m(/i 2 ) exp { - / — [1 + 7 m (#(A))] 
m(/x 2 ) 



{-//;> 

, fi^oo 0. (22) 

Here 7 m (5 f ) is a perturbative quantity, much like (3(g), which therefore vanishes as \i increases 
and g(fi) decreases. Thus, as indicated, the effective mass vanishes compared to the renor- 
malization scale when that scale diverges, and the perturbative theory becomes, effectively, a 
massless theory. 
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An important point for the phenomenology of QCD |32j is that the "light" u and d quarks 
have masses of just a few MeV at scales of /i ~ 1 GeV, so that 



m U)d (//) < A (23) 

for any scale right down to yu's characteristic of strong coupling. For any but the longest 
scales, the same holds for the s quark. Thus, for these quarks, QCD is effectively a theory of 
massless particles at any scale where we can hope to do perturbation theory. The other three, 
"heavy" quarks, c, b and t have masses whose running must be taken into account. Depending 
on the quantity being considered, each requires special treatment, and correspondingly offers 
a rich variety of theoretical and experimental challenges and opportunities. For the most part 
here, however, we shall concentrate on QCD as a theory of massless quarks, and try to address 
the question of how to use its asymptotic freedom in the "real world" of Minkowski space. 



1.4 Wick rotation 

By itself, asymptotic freedom is a striking result, and very supportive of the partonic picture of 
scaling in deeply inelastic scattering. But the problem arises of how to separate short distances, 
represented by the Born cross section in the parton model, from long, represented by the parton 
distributions. How can we increase our confidence in the heuristic arguments given above? The 
operator product expansion is designed to organize the relationship between short distances 
and long, and indeed, it will play an important role in our analysis. This analysis, however, 
now will be somewhat more complicated, just because we want to discuss cross sections in 
Minkowski space, rather than Green functions in Euclidean space. 

Since our aim is to use asymptotic freedom in Minkowski space, it is worthwhile to recall 
the method of Wick rotation, which is used to define Green functions as they approach physical 
regions in external momenta. This technique is a very general one, which applies to any 
relativistic field theory. For simplicity, consider a self-interacting scalar field. The generating 
functional for its Green functions may be defined in Euclidean space as the parth integral 



Ze[J\ = J[d4>\ exp 



d x 



) 2 + mV 



+ V((f>) + J<P 



(24) 



where V(<p) is a potential. 

Formulating the field theory in Euclidean space has a number of advantages. The path 
integral rapidly damps contributions from large values of the field. Correspondingly, the pertur- 
bative Green functions found from Ze[J] are real functions, free of singularities for real external 
momenta. In perturbation theory, their internal propagators are of the negative definite form 



A(k) 



k 2 



k 2 



(25) 



i=i 
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Such Green functions, however, cannot describe physical processes. Indeed, the S-matrix el- 
ements necessary to define cross sections are themselves found from singularities in Green 
functions when external momenta approach the mass shell, pf = m? . Wick rotation makes the 
connection between these two regions. We analytically continue Ze by replacing the real "time" 
component £4 by a complex number r in terms of two real parameters Xq and 9 as follows: 



2x4 



T 



(26) 



For arbitrary values of 9 we define a new generating functional Zq by replacing dx^ by idx^e 

/ o , \ 2 



Z g [J] = / [d<f>] exp 



-ze 



dxg d X 



1 



„2i0 



(V0) 2 + m 2 (j) 2 



+ V((f>) + J</> 



(27) 

Note that (p6|) is not a change of variables, and that Zq is not equal to Ze except for 9 = 7r/2; 
rather, we regard it as an analytic continuation of Ze in r. We must assume that such a 
continuation is possible for the full generating functional, as turns out to be the case for its 
perturbative expansion, for which the propagator is given at arbitrary 9 by 



A(k, 9) 



_p -2i0 _ P 



(2? 



This expression has no singularities for real kg, hi, so long as 9 remains finite. 



The physical field theory is generated by the limit 9 —>■ + , in which 9 approaches zero 
from above, 



Z [J] = / [d(f>] exp 



% J dx d 3 x (- ((9 O 0) 2 - (V0) 2 - m 2 2 + ietf) - 1/(0) - «V) 



In this limit A(fc, + ) becomes precisely the normal Feynman propagator, 



(29) 



(30) 



whose singularities are regulated by an u ie prescription". 

Wick rotation for a two-point function. As a practical example in perturbation theory, 
let's consider the scalar equal-mass two point function at one loop, fig. [|, directly in Minkowski 
space. Working in n dimensions, combining denominators by Feynman parameterization and 




Figure 4: Scalar two-point function. 
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completing the square in the loop momentum, we find 

d n £ 



I(P 2 



g 2 ^ 



2 2e 

9 V 



2 2e 

9 A* 



(27r)" (£ 2 - m 2 + ie)((p - t) 2 - m 2 + ze) 
d n £ r 1 dx 



(2n) n Jo (£ 2 - 2xp ■ £ + xp 2 - m 2 + ie) 2 



i f d nfJ> 

dx 



io J (27r) n (£' 2 + x(l - x)p 2 - m 2 + ie) 2 
The factor /i 2e keeps the coupling dimensionless in n dimensions. 



(31) 



In Minkowski space, with 



fi2 




V , the integrand in ( pT| ) has poles at 



±y £' + m 2 — x(l — x)p 2 =F ie . 



(32) 



The integral, however, is defined by the ie prescription, which keeps the poles off the real axis, 
so long as p 2 < 4m 2 . For these values of p 2 , the poles are in the standard arrangement shown in 
fig. [5], and the £' integration contour may be rotated to the vertical (imaginary) axis, according 
to 



pi 



where i£' n = 



e~ < 9 < tt/2 , (33) 
w / 2 is the imaginary energy along the vertical axis in the £' plane. This rotation 



/ / 

n 


/I* 

A 


X 






1 



Figure 5: Wick rotation in the energy plane. 

is in just the opposite sense from the rotation of time in eq. (p5|). For p 2 = —Q 2 < (in fact, 
for p 2 < 4m 2 ), the integral may be rotated as an analytic continuation into a purely Euclidean 
integral, 

id n £' 1 



H-Q 2 



2 2e 

9 V 



dx 



(2tt)" {1% -x(l -x)Q 2 -m 2 ) 



g 2 ^iT(2 



n, 7T 



n/2 



dx 



2 J (2n) n Jo (x(l -x)Q 2 + m 2 ) 2 -™/ 2 
where d n £' = d£' n d 3 £'. Exanding around n = 4, and defining 

e = 2 -n/2, 

4 This e is one-half of the convention used by Collins in his lectures. 
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(34) 



(35) 



we find 



Hp 2 ) 




e (4tt) 



¥ [l-e( 7E -ln47r)] 



(4tt) 2 7o 



^ 2 /■ 



i 



dx In (m 2 - x(l - x)p 2 )/fi 2 + 0(e) . (36) 



In the full theory, the pole term is removed by ultraviolet renormalization. We are interested, 
rather, in the finite, momentum-dependent remainder. It is real for p 2 < 4m 2 and is defined 
by analytic continuation for p 2 > 4m 2 , where, on the real axis, it is complex but finite, as long 
as m 2 > 0. Real, negative p 2 is equivalent to Euclidean space in which, as observed above, 
all Green functions are (relatively) real, analytic functions of their external momenta. Our 
example shows how branch points develop in their analytic continuation. Note in particular 
that the branch point occurs at "threshold" for the production of two on-shell particles. This 
correspondence of singularities to physical processes is quite general, as we shall see presently. 

Recalling that in QCD light quark masses are small compared to the natural scale of the 
coupling, we may take a special interest in the m 2 — > limit of ( j36[) . The branch point moves to 
the origin in the p 2 plane, where the function I(p 2 ) diverges logarithmically. This divergence has 
quite different interpretations in Euclidean space, where it occurs at a single point, p M = 0, and 
in Minkowski space, where it is on an entire two-dimensional subspace, the light-cone p 2 , = p 2 . 
Such divergences in Green functions will recur in QCD, and get worse at higher order. This is 
a clear difficulty for any program to exploit asymptotic freedom - since it indicates that Green 
functions have, in general, divergent sensitivity to low-momentum, long-distance physics. In 
fact, since in Minkowski space, small invariant mass does not necessarily mean small energy, the 
situation is much more complicated than in Euclidean space. It is this long-distance behavior 
whose analysis, which is still an ongoing topic of research, we shall describe in the following. 
Before we do so, however, it will be useful to intoduce a simple rule for identifying quantities 
that are not sensitive to long-distance, nonperturbative physics. 

1.5 Infrared safety 

Asymptotic freedom is useful for quantities that are dominated by the short- distance behavior 
of the theory [§]]. Such quantities, which are termed infrared safe, cannot depend sensitively 
on the masses of quarks, nor can they suffer from infrared divergences of the sort identified in 
the example above. Infrared safety is one of the fundamental concepts of perturbative QCD, 
and it makes essential use of the renormalization group. To see how, and to see further why an 
analysis of infrared divergences is necessary, we consider a generic physical quantity (say, a cross 
section) r(Q 2 / fi 2 ,a s (fi 2 ),m 2 / fi 2 ), where Q represents "large" invariants, much greater than A, 
m represents light quark masses (and the vanishing gluon mass), and fi the renormalization 
scale. We assume that r has been scaled by an overall factor of Q, to make it dimensionless. 
Now because r is physical, it cannot depend on /x, 




(37) 
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This just says that we may, if we like, expand in the coupling at the scale of the large momenta 
of the problem, and use asymptotic freedom. In general, however, we pay the price of introduc- 
ing dependence on the large ratio Q/m, which typically occurs in powers of logarithms. The 
presence of such logarithms (an example of which was illustrated above) can make the pertur- 
bative expansion unusable. On the other hand, pertubation theory can be used if r happens 
to be infrared safe. To be specific, we shall demand that r behave in the large \i limit as 

r^,a s (^,^^j^oo r(^,a s (^+O^ry a >0. (38) 

That is, r should approach a limit as m/fi — > 0, with Q/fi held fixed, with corrections that 
vanish as a power of m. Again, I{p 2 ) discussed above, is an example of an infrared safe quantity 
so long as p 2 ~ Q 2 ^ 0. 

Much of the remaining discussion will center on identifying infrared safe quantities, and 
separating them from long-distance dependence. Indeed, this is the essence of the QCD jus- 
tification of the parton model. In this connection we note that infrared safety can apply not 
only to cross sections and other direct physical observables, but also to renormalization-group 
variant quantities that obey equations of the form 

(^-Tr)r = 0, (39) 

with 7r an anomalous dimension. 



1.6 QCD in the n plane 

We'll end this introductory section with a few comments on the double role of dimensional reg- 
ularization in perturbative QCD. Dimensional continuation serves to regulate both ultraviolet 
(n < 4, e > 0) and infrared (n > 4, e < 0) divergences. When the number of dimensions is 
less than four, the volume of phase space is decreased. Because ultraviolet divergences result 
from the large number of states at high energy, they are softened as the number of dimensions 
is decreased. In contrast, infrared divergences result from singularities in the integrands of 
momentum space integrals, and when these singularities are spread out over a larger phase 
space at higher dimension, they are softened. The question arises, however, how dimensional 
continuation can handle both problems at the same time. 

To see how it works, we recall the basic path from a Lagrangian to cross sections and 
other physical quantities. This is illustrated by the sequence below: 

Cqcd - G (re9 \ Pl ,... Pn ), n<A 

- G^ n \ Pl ,...p n ), n<4 + A 

-> S {unphys \p u ...p n ) , 4<n<4+A 

-> T^ unphys \ Pl ,... Pn ), 4<n<4 + A 

-> T ( P h ^(p 1 ,...p n ), n = 4. (40) 
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The first step is to generate regularized off-shell Green functions, denoted from the Feyn- 

man rules of the theory. As we shall see, these Green functions are free of infrared divergences, 
and so are rendered completely finite for n < 4. Here is where the theory is renormalized, 
producing a set of renormalized Green functions, G^ rer ^ which are now analytic functions of 
n, treated as a complex number, in some strip in the n-plane, 4 — A' < n < 4 + A, with 
A, A' > 0. It is thus possible to analytically continue the renormalized Green functions to n 
slightly greater than four. Because the infrared divergences of the theory are regulated in this 
region of n, we can now define S-matrix elements, S^ unphys ' for the renormalized theory. These 
S-matrix elements, although fully renormalized are nevertheless unphysical, because they are 
only infrared finite for n > 4. That is, they are not infrared safe. They can, however, be used to 
compute infrared safe quantities in n > 4, T ( un P h v s ) ; which, finally, may be continued to n = 4 
to derive finite, physical predictions, T^ phys ^ from the theory. It is important to stress that QCD 
for n > 4 in not physical QCD, and cannot be used to compute physical quantities unless they 
are infrared safe, and thus have finite limits for n — > 4. Having made these observations, we 
are ready to begin our analysis of infrared divergences in field theory. 



2 Long and Short Distances in Minkowski Space 

In this section we shall discuss how to analyze and classify sources of long-distance behavior 
in perturbation theory. We do so not because we expect perturbation theory to be correct 
at long distances, but rather to identify and eventually calculate short-distance (infrared safe) 
quantities for which perturbation theory may reasonably be trusted. 



2.1 Example: IR and CO divergences in the massless vertex func- 
tion 

A very informative example, which already illustrates many of the general properties we will 
identify below, is the fully massless three-point function at one loop, with two on-shell external 
lines. We shall start, as above, with this vertex in a scalar theory, shown in fig. ||a, 

3 3e f / \ 

A = 9 V J (2<r (P + ie)((pi - k) 2 + xe) {{p 2 + k) 2 + ie) ' (41) 
This integral may be evaluated just as the two-point integral above, by Feynman parameter- 
ization and n-dimensional integration. The resulting (two) parametric integrals can then be 
done in terms of Euler beta functions, with the result, 

f 1 g 2 ( 4vru 2 V / ,B(-e,l-e) . . 

/A = HS// - JTyi 2 ■ r 1 + 6 " ' 42 

q 2 (47r) 2 \—q 2 — ie) — e 

with q 2 = 2pi -p2, and where again, e = 2 — n/2. This scalar integral is ultraviolet-finite, so the 
double pole in e is entirely infrared in origin. Before discussing it further, we may exhibit the 
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Figure 6: (a) scalar vertex and (b) gauge theory electromagnetic vertex. 

corresponding result for the electromagnetic vertex function in a massless gauge theory, found 
from the diagrams of fig. ^3. 

For zero-mass fermions (quarks) the electromagnetic vertex is given in terms of a single 
form factor. For example, in a quark-antiquark production process we have 

r M(<? 2 , e) = -iefi e u (pih IJ ,v(p 2 ) p(q 2 , e) . (43) 
At one loop the form factor is given by 

<^-^) +4 )> (44) 

where CV = (iV 2 — 1)/2N (see Appendix A). Up to this group factor, the leading, double pole 
term is essentially the same as in the scalar case, eq. d42|) . These double poles are common in 
dimensionally regulated massless integrals. To understand their origin, we may return to the 
scalar case. 

Consider the integral I a over a region of momentum space where the loop momentum k 
is small enough that k 2 <C pi ■ k, p 2 ■ k. Neglecting the k 2 terms compared to Pi ■ k is called the 
"eikonal" approximation. It is a subtle approximation in Minkowski space, where Pi ■ k small 
does not necessarily imply that k 2 is smaller. Accepting that we shall have to return to this 
point later, we work in the eikonal approximation in a frame where 

Pl = (p+,0-,0 ± ), P2 = (0,pZ,0±), (45) 

and where plus and minus components are defined by ± f 3 ), so that v 2 = 

2v + v~ — v 2 . 

The integral becomes 

7 (cik) J_ f dk + dk~d 2 k ± 

A ~ 2q 2 J (-k~ +ie)(k+ + ie)(2k+k~ - k ± 2 + ' 1 ' 

In this integral we easily identify three limiting regions which lead to logarithmic divergence. In 
the first, all four components of k^ vanish together; this we shall call the "soft" region. In the 

19 



p{q\e) 



a, 
'27r' 



C f 



47r/i 2 y r 2 (i-e)r(i- 

-q 2 -ie) T(l - 2e) 



other two, the component of k^ parallel to either p\ or p 2 remains finite, while the remaining 
components vanish in such a way that k + k~ ~ kj_; these we refer to as "collinear" regions. 
Momentum components in these regions are of the order of y/q 1 times powers of a "scaling" 
variable A, which vanishes at the points in momentum space where Ja is singular, 

1. k»~\^z "soft" 



2. fc ± ~ yfq 1 (47) 
k T ~ \ 2 \f(p "collinear" 
kl-X^q 1 . 

The lograrithmic divergences in Ja may be made explicit by changing variables in each of these 
regions to A and a set of scaled momenta k^ 1 = k^/X 0- , with a the power appearing in eq. (4~7). 

The collinear singularities of (^) are the characteristic feature of infrared sensitivity 
in Minkowski space, in which on-shell lines need not have vanishing momenta. Clearly, this 
complicates the situation relative to Euclidean space. 

While soft and collinear divergences were relatively easy to identify in the simple example 
above, it is still natural to ask how to identify infrared sensitivity at higher orders, and in other 
processes. To answer these questions requires a more general analysis, to which we now turn. 



2.2 Analytic structure and IR divergences: Landau equations and 
physical pictures 

Our aim in this section is to systematize infrared analysis in Minkowski space [ |TT| , O . In this 
way, we shall find it possible to separate differing momentum scales and to apply the operator 
product expansion in the presence of light-cone singularities. 

Let us return to the massless scalar triangle again, treating it this time from a more general 
point of view. Introducing Feynman parameters, the scalar triangle is given by (suppressing 
the coupling), 

d n k f 1 da\da 2 da 3 5(1 — Yh=i a i) 

, (48) 



(27r) n Jo D 3 
where the new denominator is 

D = aik 2 + a 2 (pi - kf + a 3 (p 2 + kf + ie . (49) 



Again, the scalar integral in eq. fl48|) is ultraviolet finite, so the poles in (^) at e = must 
come from infrared sensitivity due to the vanishing of D. 

We shall make strong use of the analytic structure of the integral Ja, which is defined in 
terms of integrals in the complex k^ and planes, and whose singularity structure is defined by 
the u ie" in eq. (f49D . Since D is quadratic in the momenta, it has no more than two poles in any 
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momentum component, when the others are held fixed. When these poles are at real values, 
D vanishes at two points along the contour. Now we recall that by Cauchy's theorem such a 
contour integral may be deformed between any pair of paths so long as this deformation crosses 
no singularity, in this case no point D = 0. Therefore, so long as the solutions to D = are 
separated, the relevant contour may be deformed to make the integrand bounded at all points, 
leading to a finite result. On the other hand, if the poles coalesce (they are automatically in 
opposite half-planes), the contour can no longer be deformed, and the result may be singular. 
This is called a "pinch" of the contour. Because D is quadratic in momenta, this is equivalent 
to the condition 

D( ai ,k^p a ) = 0, (50) 



at D = 0. Assuming that we may choose to do any of the momentum integrals first, a necessary 
condition for a singularity is to have a pinch in every loop momentum component. 

Similar considerations apply to each of the aij's. In this case, however, D is linear in 
each «j, so there are never two poles to pinch, and the a.^ contour may always be deformed 
away from a pole, except at the origin, since this is where its integral originates. A pole may, 
however, migrate to an end-point a, = 0, or it may be that at D = 0, if — mf = on the line 
corresponding to a iy so that D is indepedent of CKj 

These requirements apply equally to any diagram at any order, with line momenta ii and 
loop momenta k s , and may be summarized as 

either if = mf, or = 0, 

and a ^i e is = 0, (51) 

i e loop s 

for all i and s. Here ej S is an "incidence matrix", which takes the values +1 and —1 when 
line momenum ii flows in the same direction or opposite direction as loop momentum k s , 
respectively, and is zero otherwise. Eqs. (|5T|) are commonly known as the Landau equations 



For the three-point vertex functions of fig. || the Landau equations are simply 

aik" - a 2 (pi - k f + a 3 (p 2 + fc)" = . (52) 
One solution to these equations corresponds to the soft limit of vanishing gluon momentum, 

F = 0, (a 2 /a 1 ) = (a 3 /ai) = 0. (53) 

Another set of solutions correspond to the collinear limits, where k becomes proportional to pi 

or p 2 , 

k = (pi, a 3 = 0, ai( = a 2 {l - C) 

k = -C'p2, « 2 = 0, a x C = a 3 (l - CO . (54) 
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Are these the only solutions? The task of tracking down pinches is greatly simplified by an 
observation due to Coleman and Norton [34]]. We begin by identifying the products a^i for 
each on-shell line with a space-time vector, 



aniS = Axf . (55) 

Suppose we identify further = Ax®/£® as the ratio of the time component of Axi to the 
energy Then 

Axf = AxX, , (56) 

with 

A 1 

the four- velocity of a particle of momentum . With this interpretation, Arc, may be thought of 
as the four- vector describing the free propagation of a classical on-shell particle with momentum 
£i. In an alternate derivation of the Landau equations, based on time-ordered perturbation 
theory, this physical picture emerges automatically. This derivation is given in Appendix B. 

Let us apply this analysis to the collinear pinch of the triangle diagram, eq. (|53D , where 
k becomes collinear with pi. Consider the vectors associated with the two on-shell lines at this 
pinch, 

Aa&.a = a 2 ( Pl - hY = aiF = Ax k " . (58) 

They are equal. Now consider the diagram shown in fig. |7|, We have contracted the single 
off-shell line p 2 + k to a point, corresponding to = 0, i.e., no propagation for this line. Any 

Pi 

Pi - 
7 

Figure 7: Reduced diagram corresponding to a collinear pinch surface. 

such diagram, in which off-shell lines are contracted to points, is called a "reduced diagram". 
The reduced diagram for this pinch describes a physical process, in which two on-shell massless 
particles, of momenta k and p± — k, are created at vertex 1, and propagate freely to vertex 2, 
where they combine to form the outgoing massless particle p±. This is kinematically possible 
only because the lines are massless. We have found that this collinear pinch surface describes a 
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physical process, in which vertices may be identified with points in space-time, between which 
particles propagate on the mass shell. 

The generalization of this result to a completely arbitrary diagram is quite straightorward, 
and only requires the schematic subdiagram shown in fig. |3|, in which a loop of n lines is shown. 
For this loop to describe a portion of a physical picutre, the four-vector separation across any 




Figure 8: Schematic loop in an aribitrary reduced diagram. 

line must equal the separation derived by going around the other n — 1 lines on the loop, 

Ax 12 + Ax 23 + ... + Ax nl = . (59) 

Given the identification in eq. J55D, we see that this requirement is identical to the Landau 
equations, fl5"T|). 



2.3 Power counting and pinch surfaces 

In the Landau equations and the physical picture analysis, we have powerful tools for the iden- 
tification of sources of long-distance sensitivity in perturbation theory. But even a pinch is 
only a necessary condition for infrared divergences. In many circumstances, the perturbative 
integration contour may pass through a pinch surface without producing a singularity. The con- 
tributions from such regions may be vanishingly small in the limit of large momentum transfer, 
consistent with the requirements of infrared safety. When summing perturbation theory to high 
order, such regions may again become important, but for the present we shall look for a further 
necessary condition for infrared divergences at finite order. To find it, we shall study how to 
bound integrals near pinch surfaces. We refer to this process as "infrared power counting" , in 
analogy to the ultraviolet power counting employed in perturbative renormalization . 



A general pinch surface is depicted schematically in fig. [| At each point on a pinch 
surface S, we identify coordinates that lie in the surface, which we refer to as "intrinsic", and 
those that parameterize directions out of the surface, which we call "normal" . By construction, 
the integrand is a singular function of the normal coordinates only. 

It is possible to bound the integral near S using the power counting technique illustrated 
in Sec. [2.1| for the triangle diagram, fig. || We put bounds on the ratio of the volume of normal 
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Figure 9: Schematic pinch surface S with one normal and two intrinsic coordinates. 



coordinates to the magnitude of the integrand at the pinch surface. In general, the larger the 
volume of the normal space, the less singular the integral. This power counting will have a 
dual purpose. First, it is used to put bound on integrals, and hence to identify infrared safe 
quantities, and by the same token it may also be used to identify regions in momentum space 
that may give rise to infrared divergences. 



i. We redefine each of the normal variables Kj in terms of a scaling variable A according 
to 

Kj = \ a i K '.. (60) 

We will determine the behavior of the integral when A vanishes for fixed values of the 
ratios Kj/Kj>- For instance, in the triangle diagram, eq. (|4"ID, the scalings of eq. ([47]) 
specify that the the propagators (pi — k) 2 and k 2 are linear in A, while (p2 + k) 2 is zeroth 
order in A in the collinear region where k is in the pi direction. Similarly k 2 is quadratic, 
and both (pi — k) 2 and (p2 + k) 2 are linear, for the soft scaling of eq. fl4"7|). 

ii. Given a set of powers aj, we retain only terms of lowest power \ Ai in A for each 
perturbative denominator k 2 (nj, A) — m 2 , 

k 2 {K v \)-m 2 = \ A ^f{K' ] ) + .... (61) 

We call the resulting integral over the normal variables Kj the homogeneous integral. If 
the homogeneous integral is independent of any normal variable kj, its scaling power aj 
may be reduced until that variable appears in the homogeneous integral [FT]. 

iii. The homogeneous integral for pinch surface 5* is proportional to A ns , with ns given 
by 

n S = J2 a i ~ J2 A i + S l ' ( 62 ) 

i i 

where s/ represents the (possible) power of A from momentum factors in the integrand 
(these will be important for QCD). If ns > 0, the integral is finite when all of the normal 
variables vanish according to the scaling (^1|). If ns = 0, the integrand may diverge 
logarithmically, while if ns < 0, it may diverge as a power. 
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iv. Finally, check for pinch surfaces in the homogeneous integral. If the only pinch 
surface is the original one, at which all normal variables vanish, the bound is complete 
for the singular surface S. Should the homogeneous integral have further pinch surfaces, 
however, bounds must be found for these special regions as well. These will correspond 
to subsurfaces of S where subsets of the normal variables vanish faster than others Jll| . 



2.4 Pinch surfaces for the all-order EM form factor 



We can now apply the Landau equations and power counting to identify the sources of infrared 
divergence in the electromagnetic form factor to all orders. Consider pair creation, in which 
a vector particle of momentum q, q 2 > (Z or virtual photon), decays into a massless quark 
and antiquark with no radiation. The possible reduced diagrams associated with pinch surfaces 



are remarkably simple, and are all of the form shown in fig. 10. The general reduced diagram 




Figure 10: General reduced diagram for pair creation. 



corresponds to a set of physical processes in which the decay of the vector is followed by the 
formation of two "jets" , labelled J and J' in the figure, of virtual particles in the same direction, 
and with the same total momenta, p and p', as the two final state particles. The only interaction 
between the two jets is via zero-momentum "soft" particles, labelled S. Higher-order off-shell, 
short-distance, contributions reside in subdiagram H, adjacent to the decay vertex. No other 
physical processes, in particular no finite momentum transfers between the two jets, can be 
realized, because once the jets are formed at the point represented by H, they travel apart at 
the speed of light and can never meet again at a point in space-time. 

Power counting for the pinch surfaces corresponding to the reduced diagrams of fig. |10 
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is straightforward. A natural choice of normal variables is (i) all four components of each 
loop momentum internal to the soft subdiagram S, and of each loop momentum that links the 
soft subdiagram with the jets; (ii) the total invariant mass £ 2 and the scalar product i ■ p for 
each loop I in J, and similarly for J'. All of these variables may be scaled with Oj = 1 in eq. 
fl60P . The resulting homogeneous integral has no pinch surfaces, except those that correspond 
to reduced diagrams of the same form as fig. |l(] Jl|, [TTJ] . The scaling powers of lines are Aj = 1 
for lines internal to J and J', and Aj = 2 for soft lines in S. Self-energies on jet lines need not 
be considered explicitly, because each full two-point function G 2 (k) has only a single pole at 
k 2 = 0. 

Consider an arbitrary pinch surface T with Nj lines and Lj loops in J, and similarly for 
J' and S. We assume in addition that Mj soft lines attach 5* to J, and Mj, soft lines attach it 
to J'. The power counting measure is then, according to eq. fl62|) , 

n T = A(L S + Mj + Mj, - 1) + 2Lj + 2Lj, - 2N S - Nj - Nj, + s T , (63) 

with St the numerator suppression factor. The contribution of all soft lines, loops and numer- 
ator momenta internal to 5* is just the dimension of S. For technical simplicty, we shall assume 
that the jets attach to H each by a single line, and that the soft lines are all gluons, attached 
to the jets only at three-point vertices. The dimension of S, including its external lines, is then 
4-3M/-3M//, and 



n T = Mj + Mj, + 2Lj + 2Lj, - Nj - Nj, + s' T , (64) 

where s' T meaures the power associated with numerator factors that come from jet lines and 
vertices. 

It is not difficult to find a lower bound for the numerator suppression factor s' T . We note 
that there is a factor of numerator momentum for each three-point vertex in the jets, and that 
each of these momentum factors will have to combine to form an invariant. Each such invariant 
will scale as A if it is the scalar product of two momenta within a jet, but as A if it involves 
momenta from different jets. On the other hand, the only way the latter may occur is if each 
factor of momentum is contracted with the spin tensor of one of the soft gluons. We conclude 
that 

s' T > \{V^j + V 3 ,j, -Mj-Mj,), (65) 



where V^,j is the number of three-point vertices internal to jet J. Substituting (|65j) into (p4|) 
we have 

n T >\{Mj + Mj,) + 2Lj + 2Lj, - Nj - Nj, + ^{V 3)J + V 3)J ,) . (66) 
Now for each jet we have the Euler identity 

Lj = Nj-Y,V ly j + l, (67) 

i 
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where Vi is the number of i-point vertices in the jet (counting the hard part if as a one-point 
vertex for each jet), as well as the relation 

2Nj + M J +\ = Y t iV i , (68) 

i 

which takes into account the two external lines in each jet, one attached to H, and one in the 
final state. It is now a rather straightforward exercise to show that 

n T > , (69) 



which shows that at pinch surfaces like fig. [10], the integral is at worst logarithmically divergent. 

In addition, a closer look, using the same reasoning, shows that ny = 0, corresponding to 
logarithmic divergence, can occur only when the reduced diagrams for each jet are themselves of 
the forms of standard perturbative diagrams, with only three- and four-point vertices internal 
to the jets. A five- or higher-point vertex in a reduced diagram always leads to suppression in 
the infrared. In addition, we also find that S is connected to the jets by soft gluons only; the 
interaction of soft quark lines with jet lines is similarly suppressed, although soft quark loops 
are possible internal to S. It is clear that attachments of soft lines of any flavor to H will be 
suppressed relative to diagrams where they are attached only to jets. 

In the above, we assumed that the jets attach to if by a single line. This assumption 
depends on the choice of gauge. In an axial gauge n ■ A = 0, the form of the propagator 

°^ n) = [-"i" + ^k + ^k - -J^W) WTVe ' (70) 

leads to a suppression whenever it is contracted with the momentum of the gluon itself. That 
is, the combination 

k ■ n (k ■ n) 

has no pole at k 2 = 0, and hence does vanish except at k — 0. This leads to a contribution to 
s' T of 1/2 in (|6~5f) whenever k is a jet line, even if the line k attaches a jet to the hard part H. 
Then it is easy to verify that in an axial (or other physical) gauge, reduced diagrams that give 
infrared divergences have only a single line connecting each jet to H, as assumed above. 

In covariant gauges, such as the Feynman gauge, there is no such suppression, and multiple 
collinear gluons may attach J or J' to H. This is not quite the complication it appears to be, 
however, for the following reason. At the pinch surface, we can have tlt = only when each of 
these gluons appears in the combination 

k^(k) = ^, (72) 

in which the gluon propagator is replaced by an unphysical "polarization" vector k u , equal to 
the gluon's momentum. Such gluons are sometimes described as "scalar-polarized". As we 
shall see, scalar-polarized gluons decouple from the hard scattering |J. 
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FGWM) = " TiT^ > ( 71 ) 



In summary, reduced diagrams associated with logarithmic divergence in the electromag- 
netic (or other electroweak) vertex are characterized by a simple two-jet structure, with a 
single hard-scattering function, in which the jets are connected to the hard part by single lines 
in physical gauges, and in which the soft part is connected to the jets by soft gluons only, and 
not to the hard part. 

Most of these observations are much more general than the electromagnetic form factor 
from which they have been derived. In the following, we shall apply the tools we have developed 
to a wide set of processes. 



3 Short Distance Cross Sections and Unitarity 

We are now ready to discuss a class of physical quantities that can be proved infrared safe by 
a combined analysis of their analytic structure and power counting. They are primarily cross 
sections initiated by timelike electroweak currents, the prime examples being the total and jet 
cross sections in e + e~ annihilation |3J. Other examples include the decay width of the Z and W, 
and various event shapes defined to describe cross sections and decay rates of this type. The 
analysis of these cross sections is related to the operator product expansion, but gives many 
results that cannot be derived directly from the operator product expansion. 

3.1 Cut diagrams and generalized unitarity 



To discuss cross sections of this type, we introduce the "cut diagram" notation shown on the 



left-hand side of fig. 11. This diagram represents the amplitude for a process in which a set 





Figure 11: Cut diagram and unitarity. 



of particles with momenta Pi, ■ ■ ■ p n scatter into a set £i, . . . £ m , times the complex conjugate 
amplitude for the latter to scatter into a third set k±, . . . k n r, integrated over part or all of the 
phase space of the intermediate set. Perturbation theory rules to the left of the cut are the 
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usual ones, those to the right their complex conjugates. Each such cut, which we label C in 
fig. [11], specifies a distinct intermediate state. We shall denote a particular cut diagram found 
from uncut diagram G by cut C as Gc- 

Fig. [TT| as a whole states a very useful theorem satisfied by cut diagrams. The sum over 
all cuts with fixed external momenta pi, . . . p n and k\, . . . k n > is given by twice the imaginary 
part of (— iG), the uncut diagram, 

]T G c (Pi, kj) = 2 Im {-iG{pi, kj)) . (73) 

all C 

This result is a generalization of unitarity expressed in terms of the T matrix, T = — i(S — 1), 
with 5* the scattering matrix, 

TT f = -%{T - T f ) . (74) 

It is, however, more general, because it applies for fixed spatial momenta, in both the external 
and internal loops of the overall diagram G. This result is proved in Appendix B by use of 
time-ordered perturbation theory. 

3.2 Infrared safety for inclusive annihilation and decay 

The simplest physical application of the analyticity/power counting analysis is to the total 
cross section for e + e~ annihilation or Z decay, fig. [12]. By the optical theorem, a special case of 




Figure 12: Unitarity applied to the total annihilation cross section. 

generalized unitarity, eq. ([73]), the total cross section is proportional to the vacuum expectation 
value of the time-ordered product of two currents, 

^ + t(Q 2 ) = - 2 lm^), (75) 

where the function tc is defined in terms of the two-point correlation function of the relevant 
electroweak currents (with their couplings included) as 

<q 2 ){q»qv - q 2 g^) = i f d A xe^ x < o| t j^x)j u (o) |o > . (76) 
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The proof that 7r(g 2 ), and hence its imaginary part, and therefore the total cross section, 
is infrared safe requires only that we recognize that 7r represents a forward scattering process. 
There are no physical processes in which an off-shell photon (or on-shell Z) can decay into a set of 
on-shell particles that propagate freely and then annihilate to form another photon of the same 
invariant mass. The set of particles originating from a point will receed in different directions, 
and can never meet again by physical propagation. This eliminates pinch surfaces with finite 
momentum particles, and hence any infrared divergences associated with decay processes. This 
is a particular example of the famous Kinoshita-Lee-Nauenberg (KLN) theorem . 

It is still possible to have pinch surfaces involving only massless particles, coupled to a 
single hard-scattering, as in fig. |H| Dimensional counting shows that such pinch surfaces are 
finite order-by-order in perturbation theory. They will, however, lead to an important connec- 




Figure 13: Pinch surface with only zero-momentum lines (subdiagram S) coupled to a point-like 
hard part H. 

tion of perturbation theory, even without explicit masses, to the operator product expansion, 
as we shall see in Section 7 below. 



3.3 Jet and weighted cross sections 

The class of infrared safe cross sections is by no means exhausted by the totally inclusive pro- 
cesses above. More detailed information on final states is available in jet and other "weighted" 
cross sections. 

Jet cross sections measure the probability of producing states that are identified as jet- 
like, according to infrared safe criterea. The simplest of these [|J is illustrated by fig. 14, which 



are defined by the energies flowing into each of two cones, 5, back-to-back along a fixed axis, 
at angle 9 from the beam direction. We may define a two-jet cross section in e + e~ annihilation 
as one for which 

^>l- £ , (77) 
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Figure 14: Cones for a two-jet cross section. 



that is, for which all the energy, up to a fraction e, is emitted into the two back-to-back cones. 
If the cones are sufficiently small, the event will "look like" two jets of nearly collinear particles. 
Such a cross section is infrared safe, and the proof is a variation of the proof of infrared finiteness 
of the total annihilation cross section. 

Consider for simplicity a two-jet cross section. The relevant pinch surfaces of cut diagrams 
are shown in fig. [T^. Each cut C of the diagram has a final state that contributes to the jet 



cross section, with one jet in each of the cones of fig. [14], for instance. The pinch surfaces 
associated with long-distance behavior are easily verified to be of exactly the same form as 



those for quark- ant iquark production, fig. [10], except that now there may be any number of 
particles from each jet in the final state. At the pinch surface, if cut C of reduced diagram 
R contributes to the two-jet cross section, then every cut of R contributes to the same cross 
section. The sum over cuts, however, may then be carried out using the generalized form of 



unitarity, fig. [IT], to derive an integral which has no pinch surface at all corresponding to the 
two-jet configuration, by exactly the same reasoning as for the total cross secction. 

The difference between total and jet cross sections is that momentum integrals for jet 
cross sections encounter boundaries between two- and three-jet events (for instance), which 
are absent in the total cross section. At such points, the phase space is discontinuous, and 
integrals may not be deformed after the sum over cuts. Manifolds of such points are of reduced 
dimension, however, compared to the general pinch surfaces of cut diagrams, which produce at 
worst logarithmic divergences [|13|, [IIJ. This reduced dimension weakens their power-counting, 
and they do not give rise to infrared singularities order-by-order in perturbation theory. Similar 
considerations apply to pinch surfaces with more than two jets. 

In weighted cross sections, final states are weighted according to "shape variables", 
S n (pi . . -p n ), which are functions of the momenta of particles in the final state. The shape 
variables may or may not be chosen to enhance jet-like configurations. A general cross section 
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Figure 15: Cut reduced diagram for two-jet configuration, 
at fixed shape variable S is of the form 

°S = J dT n Jj- SnijPl ■■■Pn), (7* 



where the S n are functions of the momenta in an n-particle final-state, whose phase space is 
denoted by dr n . A weighted cross section is infrared safe whenever the same mechanism that 
produces infrared safety for jet cross sections applies to them as well. That is, it is infrared safe 



if, whenever a pinch surface of a cut reduced diagram like fig. [15] contributes to the weighted 
cross section, the weight is the same for every possible cut of the reduced diagram at the 
pinch surface. If this holds, the divergences of individual cuts cancel, since the sum over cuts 
produces an integral without the corresponding pinch surface. This condition will be satisfied 
so long as the weight function does not distinguish between states in which one set of collinear 
particles is substituted for another set with the same total momentum, or when zero-momentum 
particles are absorbed or emitted. Quantitatively, these requirements may be summarized by 
the conditions 

S n (px . . .pi . . . p n _ x , Xpi) = <S n _!(pi ...pi + Xpi... p n _ x ) . (79) 

Perhaps the best known weight is the "thrust", defined |37| for an event with n particles of 
momenta p*j by 

1 " 

T = - max ft 22 | pi ■ h | , (80) 

| Pi I j = i 

with the maximum taken over all unit vectors h. The direction of h that produces the maximum 
is known as the "thrust axis" . 

Other important weighted cross sections iteratively assign particles into jets of momenta 
Pj. The algorithm begins with each particle defined as a separate jet. At each stage in the 
iterative process, a set of variables y^, 

y jk = -J{P j ,P k ) (81) 
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are computed for each pair of jets. The function / is chosen to be consistent with infrared 



safety. An influential choice |38| defines the "&t algorithm", based on 

y jk = min (E*, E 2 k )(l - cos 9 jk ) , (82) 

where Ej and E k are the energies of jets j and k in the overall center-of-mass frame, and 9j k is 
the angle between them. For this, and other such variables, if one or more of the are smaller 
than a fixed value y cu t, the pair whose value is minimum is combined to form a single jet and 
the process is repeated. The number of jets left when all yjk are larger than y cnt defines the 
number of jets in the event, which may be plotted against y cut and compared to perturbative 
or event-generator predictions. 



4 Factorization and Evolution in DIS 

4.1 Venturing out on the light cone 

Although the set of infrared safe cross sections is large, it is a small subset of interesting 
scattering experiments. In particular, QCD cross sections in which one or more of the incoming 
particles is a hadron are never fully infrared safe. For DIS, this is already clear in the parton 
model, in which parton distributions must be taken from experiment. This will remain the case 
in field theory, where we shall discover how to extract infrared safe quantities from many (not 
all!) cross sections, a process known as factorization 0. 

Factorization, as mentioned at the outset, is a realization of the separation of long- and 
short-distance dynamics in the presence of massless particles in Minkowski space. In this 
context, we shall derive a generalization of the parton model in QCD, and a field-theoretic 
definition for parton distributions |M ED|. Let us begin with dimensionally regularized deeply 



inelastic scattering of a light-like parton. Now our initial state is fixed on-shell, and is sure to 
include singularities in Feynman integrals assoicated with the light cone. 

Amplitudes for DIS are full of collinear singularities; to analyze them we turn again to 
a pinch surface analysis via physical pictures. That is, we look for the most general physical 
process in which a single incoming parton absorbs a spacelike photon to form a final state. 



The result, analogous to fig. [15] for e + e~ annihilation, is shown in fig. [16]. In this process, 
the incoming particle of momentum p produces a jet of collinear lines J(£), one of which, a 
parton of flavor % with momentum £p, absorbs the virtual photon (or other electroweak boson) 
of momentum q at H, producing a set of outgoing jets. All jets, including the remains of 
the incoming jet of lines collinear to p, propagate into the final state. The jets may interact 
mutually only through the exchange of soft particles. Because all particles in the physical 
picture must have positive energy, we find that 

1>£>x. (83) 
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Figure 16: Pinch surface for DIS. 



Thus, the struck particle is a parton of the incoming particle, with a proper fractional mo- 
mentum. That there is only one struck parton, at least in a physical gauge, follows the power 
counting of Sec. 2.4, especially the discussion of eq. (|71|). An exception is that in covariant 
gauges the physical parton may share its momentum with unphysical, scalar-polarized gluons. 
We shall incorporate this possibility below. 

Next, if we sum over all final states at fixed q, soft and collinear divergences associated 
with the final-state jets and their interactions cancel, by the same unitarity arguments as for 
e + e~ annihilation. This may be verified by using the relation of the hadronic tensor to 
the forward Compton amplitude T^, 

W^ = 2 1mT^ (84) 

with 

T,u = ^~ I d 4 xe^ < h(p)\ T J^JyiO) \h(p) > . (85) 

o7T J 

The absence of pinch surfaces involving final-state jets in the forward-scattering amplitude 
follows easily from the much simpler physical pictures for T^, which involve the incoming jet 
only, as in fig. [17|. Note that, as in e + e~, this result, and the cancellation of infrared sensitivity 
to the final state, applies to any weighted DIS cross section consistent with eq. flT9|), and not 
only to fully inclusive DIS. Infrared sensitivity from the incoming jet remains in all these cases, 
of course. 
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Figure 17: Pinch surface for forward Compton scattering. 

In summary, in the neighborhood of any of its pinch surfaces, the DIS hadronic tensor 
takes a form that is very suggestive of the parton model, 

= £ I' di H% a (q, &)J i>aP (0 . (86) 
i Jx 

The strategy of factorization is to separate the long-distance contributions from the hard- 
scattering part. Unlike the operator product expansion in Euclidean space, however, such a 
factorization cannot be expressed as a finite sum of coefficient functions times matrix elements, 
even to the leading power in Q 2 . This is because the short-distance part and the incoming jet 
remain in covolution form, tied together by the fractional momentum £ of the struck parton. 
The finite sum, therefore, is replaced by a convolution of functions. 

As in the operator product expansion, however, we can express the long-distance function 
in terms of the matrix elements of partons that connect the incoming hadron to the hard 
scattering. Thus, the long-distance contributions of all diagrams with pinch surfaces like those 
in fig. [Tj| when the scattered parton i is a quark are summarized by matrix elements of the 
form 

•W = JeT ^-e~^ + y- < h(p, a) | q a (0 + , y~, ± )q p (0) | h(p, a) > , (87) 

spin 
U 

with a and j3 Dirac indices. Instead of a local composite operator, we find a product of operators 
separated by a light-like distance y~ . In this fashion, the plus momentum of the parton is 
fixed, although its remaining components, to which H is insensitive, are integrated over. As 
we shall see shortly, the matrix elements of such composite operators are ultraviolet divergent 
in perturbation theory, and require renormalization. This is similar to the local composite 
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operators of the OPE. The effective ultraviolet cutoff for the matrix element is referred to as 
the "factorization scale", separating the short- distance from the jet functions. 

A Fierz projection between the hard and jet functions, followed by some power counting, 
shows that the leading Dirac structure may be projected out by the trace of J a p with a single 
Dirac matrix, for unpolarized scattering. This results in the following definition for the quark 
distribution |IU . 

1 roc dy - I 



(£, M 2 ) = ^E r -^e-^ + y < h(p, a) | q(0+, y~, ± ) -n ■ 7 q(0) \ h(p, a) > , (88) 

I „ J-oo 2,71 I 



where n M is the lightlike vector directed oppositely to the incoming momentum p, n M = 5^. 
Similarly, when the scattered parton is a gluon, we encounter the distribution 

<f>G/h(L /i 2 ) = J dy-e- 1 ^- E E < h(p, I y" 0)F"+(0) | h(p, a) > , (89) 

now a matrix element of the field strengths F^. These matrix elements are illustrated in fig. 



18a. 



In summary, the inclusive DIS hadronic tensor, and hence its structure functions, may 
be written in the following factorized form, which generalizes the OPE and which is clearly a 
justification of the use of the parton model in this process, 

'x Q 2 x 



F^ h \x,Q 2 ) = E !\0P( 

i=fJ,G Jx 5 ^ / 



(90) 



Compared to the parton model, the coefficient functions are now expansions in the strong 
coupling, and the parton distributions have become functions of the factorization scale. 



4.2 Gauge invariant distributions; schemes 



The distributions, ( p8[) for the quark, and (p9|) for the gluon, have been justified so far only 
in physical gauges, where only a single, physical parton connects the hard scattering and the 
p-jet on either side of the cut in fig. [T(| At the same time, these distributions are not gauge 
invariant. Their gauge variations, however may be absorbed in the hard scattering functions. 
This is because a change in the gauge of the matrix elements (|88|) or (j89"D is equivalent to a 
phase acting on the quark or gluon operators only. 

Alternately, we may define gauge-invariant distributions by connecting the physical quark 
and gluon fields in the matrix elements with ordered exponentials (Wilson lines) along the n M - 



light cone between the fields |]39fl , 

q(y~) n ■ j q(0) -> q(y")Pexp 



v- 

-ig I dX n ■ A(\n 
lo 



n ■ 7 q(Q) 



36 



1 




f 

J 
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(c) 

Figure 18: (a) Gauge- variant parton distribution (flavor i). (b) Gauge-invariant parton dis- 
tribution with ordered exponential in the n directions, (c) Graphical representation of the 
"eikonal" propagator and vertex. 



Y, F + ,(y-) i^ + (0) 



ig d\n- A(\n fl ) 



(91) 



where n ■ A in the quark distribution is in the fundamental (quark) representation, and n ■ A in 
the gluon distribution is in adjoint representation. The gauge invariant distributions reduce to 
fl88|) and ( j39|) in n ■ A = (often identified as A + = 0) gauge. Again, any of these definitions 
require renormalization. They are commonly defined by MS prescriptions, and are referred to 
as MS distributions. 

Gauge-invariant parton distributions have perturbative expansions, defined by the dia- 



grams of fig. p~8|b, which include propagators and vertices generated by the ordered exponentials 
of eq. (|9~ll). The combination of a propagator and a vertex, illustrated by fig. [l8|c, is given by 
the expression 

(92) 



n ■ k + ie 



-ign» (Tf ) 
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where R = F for the quark distribution and A for the gluon. The linear propagator is sometimes 
referred to as an "eikonal line". 

The MS matrix elements just defined are not the only possible choices for the distributions. 
Other choices may be defined by convolution with any IR safe distribution D, as 

which preserves the factorization (p0|). Of particular interest is a choice that absorbs the full 
MS F 2 DIS coefficient function C 2 . Defining C { 2 jMS \z) = E f Q 2 f c fj (z), we take 

DfM/v) = (v/0(c)fM/v)- (94) 
In this "DIS scheme" , the parton model relation 

J2Q 2 f x<j ) f / i :\x,^) = F 2 ^(x,^) J (95) 
/ 

in eq. (|T3|), holds by construction to all orders in perturbation theory. Note, however, that the 
Callan-Gross relation in (|i~3l) , which involves F\ as well as F 2 , cannot be exact at the same 
time, and so inherits corrections even in DIS scheme. 



4.3 One-loop distributions and coefficient functions 

We are now in a position to compute corrections to parton model relations such as (|H|) in 
QCD. The complexity of these calculations increases precipitously with order, but the pattern 
is well illustrated by one-loop considerations. 

The first goal is to compute infrared safe coefficient functions in a given factorization 
scheme, by comparing perturbative expressions for structure functions FfJ> of partons of flavor 
/ with the distributions <pi/f of parton i in parton /, using the factorized expressions eq. (pOf) , 
expanded to the appropriate order in a s . Neither F^> nor (pi/f is infrared safe, and must 
be defined through dimensional (or some other) regularization, but the resulting coefficient 
functions are infrared safe. Once the coefficient functions are determined to some order, the 
same factorized expressions, fl90D , may be applied to experimental measurements of the F^ 
to determine the physical values of parton distributions in hadrons h. This process is similar 
to the parton model, and as in the parton model, many predictions result by using the parton 
distributions in different processes. 

The one-loop corrections to quark DIS are given by the cuts of the graphs in fig. |l9|a. The 
one-loop corrections of the perturbative distribution <pf/f of a quark in a quark (flavor /) are 



given in fig. |19l b for A + = gauge. A short calculation shows that the gluon emission diagram 
in 



l£b is given in n dimensions by 
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(a) 





(b) 

Figure 19: (a) Uncut corrections to quark DIS. (b) One-loop corrections to distribution of a 
quark in a quark in A + = gauge. 



The phase space of the gluon is reduced to a transverse momentum integral by fixing the 
longitudinal momentum carried by the upper quark line to £p. The remaining integral is, as 
expected, divergent at = 0, corresponding to a collinear divergence from gluon emission. As 
anticipated above, it is also ultraviolet divergent. Such a scaleless integral must be defined to 
vanish in dimensional regularization. This does not mean that the distribution is zero, however. 
Rather, the scaleless kx integral is replaced by a ^-dependent ultraviolet counterterm, which 
removes the unphysical contribution of transverse momenta much larger than the factorization 



scale ||39|| . In this manner, we define the MS distribution at one loop as a "pure counterterm". 
Adding the contribution of the virtual diagram, which is proportional to 5(1 — £), we find the 
explicit expression 



7/. 

,(1) ; 



2tt 



+ 7*-ln47r)P$(0 



where the distribution PV is given by 



(97) 



(9* 



For reasons that will become clear shortly, pty is referred to as an "evolution kernel". The 
"plus distribution" is defined through its integrals with smooth functions /(£) by 




(99) 



Other distributions <j>i/j(Ji) of parton i in parton j at one loop have exactly the form (]97|) 



a. 



2tt 



(e, e) = £ ( — + 1e ~ In 4vr ) P« (0 + . . . • 



(100) 
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Each P-jj is one of the set of evolution kernels, 



P$l(0 = C F [2D(t)-l-t + l6(l-z)] 

pSUt) = T F [(i-o 2 +e\ 



P$ G (fl = 2C A [£D(Z) + ( ~ + e) (1 - + I - < )] - - ( 1 - < ) ■ ( 101 ) 



1 

3 



where £>(£) = [1/(1 - + . 

At zeroth order in the strong coupling, we have 

$//(*) = */'/*(!- 



(102) 



for any flavors / and /', which simply states that without interaction parton / remains itself. 
Then, at lowest order in (|90"D we recover the partonic results 



x) 



(103) 



Finally, the perturbative coefficient functions may be determined from the one-loop ex- 
pansion of eq. (pOD . For instance, at one loop, we have (using (|102|) and (|103|) ), 



F 2 {M (x,Q 2 ) = C 2 {M (x, %,a s (^)) . (104) 



At one loop, the electromagnetic structure functions of fig. [TTS| a are (for details of the calculation, 
see H 0) 



Q f 2 x 



a s 1 D (i) A . . 4tt^ 2 



-C F {{l-x 



2^ ( H 1 -X) 
\ 1 - X 



2(1 -x 



] + t)*(i-*)-(i+»") 1 _ x 



9 . lnx 

+ 3 + 2x} 



a-. 



= ^^-Q/C^a;. (105) 
The determination of one-loop quark coefficient functions in the MS and DIS schemes is now 



a matter of subtraction, using (j97|) and ( [L05|) in (|104|) for MS, and then (P3)-(|95|) for the DIS 
scheme. 



4.4 Evolution 



The coefficient functions C a in eq. ([50) depend, beyond lowest order, on the momentum transfer 
Q 2 . This "scale breaking" is a refinement on the scaling behavior of the parton model. The 
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dependence of the coefficient functions, and hence of the structure functions themselves, on 
momentum transfer can be computed. 

The factorization formulas eq. (|90|) depend, as we have seen, on a factorization scale 
[A, at which short and long distances are separated, /i may be interpreted, in turn, as the 
renormalization scale associated with the MS renormalized parton densities, matrix elements 
such as eq. (|^). In a manner analogous to composite operators in the operator product 
expansion, the parton densities have calculable /i-dependence. This remarkable result is referred 
to as "evolution". Because the coefficient functions depend on Q only through the ratio Q/[i, 
their evolution in /i determines their dependence on Q. This enables us to relate structure 
functions and parton densities measured at one scale to other scales, greatly extending the 
applicability of perturbative analysis. 

It is simplest to illustrate evolution with a "nonsinglet" distribution. An example is the 
difference between proton and neutron structure functions, 

Fa (NS) = p^) _ K (n) ? (1Q6) 

which combines the contributions of the two nucleon states weighted by (twice) their isospin. 
In this combination, the contributions of gluons and "sea" quark pairs produced in virtual 
processes cancel, leaving factorized expressions for the F a ^ NS ' in terms of "valence" quark 
distributions (val ' ) . In this case, we have the difference between p and n valence distributions 
of each flavor, but for massless quarks the coefficient functions are all the same up to factors 
of Q 2 , and we may suppress partonic indices, 

F^ NS) = [ |C X ^) (x/e,Q 2 /^ 2 ,« s (/, 2 )) <^(£,/, 2 ) . (107) 

Under moments with respect to x, defined as 

f(n)= f rfxr'/W, (108) 
Jo 

the nonsinglet structure function factors into a product of moments, one for the parton distri- 
bution and one for the coefficient function, 

A (JW) (n, Q 2 ) = C<M (n, Q 2 //i 2 , a s {^)) £M>(n, a s {^)) . (109) 

The essential ingredient in evolution is the independence of the physical structure functions 
from the factorization scale /i, 

^F 1 (^)(n,Q 2 )=0. (110) 
an 

Given (|107|) and ( |110| ), (ft and C obey the joint evolution equations, 

(ij-hit (n,a.(jA*j) =- ln (a s ( l i 2 )) = -fi^lnC^ (n,Q 2 /^ 2 ,a s (fi 2 )) . (Ill) 
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The anomalous dimension 7„ is a function of a s only, since this is the only variable which (p 
and C have in common. 

Logarithmic Q 2 dependence in moments of F may conveniently be computed from ( [111] ) 
by setting \x = Q, and solving for the \i dependence of 4>{n). Contributions from C{n) will then 
be a simple expansion in a s (Q), since all logarithms of Qj\x will vanish. The relevant solution 
for is 

^\n^ 2 ) = ^(n,/io a )e^|-ijr ta ^ /Wa dt7»(a.0ioV))J 

= *-W)«p{-£-j£ < + m W v) + -}- < 112 > 

where 7„(a s ) = (a s /7r)7^' ) + . . .. Using the asymptotically free running coupling eq. in 
the second line of ( |1 1 2| ) , we find that the resulting \x (and hence Q) dependence is logarithmic 
in QCD, which accounts for the mild nature of scale breaking. For a frozen coupling, or one 
that runs to a finite value, scale breaking is power-law in Q, in apparent contradiction to the 
successes of the parton model. This is a fundamental success of quantum chromodynamics, 
which played a central role in its acceptance. Indeed, the results described here meet the 
challenge set in Section 2 above, to use asymptotic freedom to account for the approximate 
scaling of DIS in the presence of multiple scales. 

Although the analysis we have just described is particularly simple for moments, in most 
practical cases, it is best to work with the parton distributions themselves. The full set of 



moment evolution equations (|1 1 1| ) are very conveniently summarized by the celebrated DGLAP 
(Dokshitzer-Gribov-Lipatov-Altarelli-Parisi) |7j integro-differential equation for the <pi/h(x, fx), 

dfi , u.<; J r s £ 

where the distributions summarize a matrix of "singlet" anomalous dimensions as their 
moments, 

\tr x P ii {fi,a.)=-li j (n). (114) 

The Pij(x,a s ) are power series in the strong coupling given to one loop by the distributions 
above, for instance, 



a 



P qq (x,a s ) = ^f(x) + ... . (115) 



7T 
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Eq. ( p. 13) is one of the most useful tools in perturbative QCD and in the search for new 
physics, since it enables us to connect experiments at widely differring momentum transfers, 
and to predict the outcomes of experiments even at very high energy. 



4.5 The light-cone expansion 

Before generalizing factorization beyond DIS, it is useful to acquire some insight into the field- 
theoretic content of evolution. Consider the moments of a nucleon quark distribution, taken 
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for simplicity in A + = gaugeQ (and supressing the spin average), 

- 1 /"l f°° - _L 

^a/N(n,a s (fx 2 )) = — dxx n ~ l dy-e- iy ~ xp < N(p) \ q(y-)l + q($) \ N(p) > 



2tt Jo 

1 00 1 r n 

o- E -r < ^(p) 1 \(d + rm] 7 + ?(o) 1 iv(p) > 

Z7r m=0 m ' 

1 /-oo 



(p+)n 



x / fa"" 1 / dy-(y-) m e~ iy xp 

JO '-oo 

< JV(p) I [(-i9 + )"" 1 g(0)l 7 + g(0)|iV(p) > . (116) 



In the second equality, we have formally expanded the q field about y~ = 0, that is, on the 
light cone, and in the third we have done, first, the resulting y~ integrals to get delta functions, 
and then the x integrals (treating j Q dx5(x) = 1). We see that only a single term contributes 
to the sum, corresponding to the local product of the quark field with the n — 1st derivative 
of its conjugate. Thus, moments of parton distributions are related to local operators, with 
dimensions that increase with the moment variable n. 

This rough discussion can be carried out in an arbitrary gauge, and the full set of relevant 
gauge-invariant operators found in this manner is 

/n-l \ 



O^i-Mn = £(0) ( JJ iD^[A] ) 7 ^g(0) (117) 

for quarks and 



vi=l 



/n-l 



C Mi-Mn = F^ a (0) Y[ iD^[A] F Q ""(0) (118) 



\i=2 



for gluons, with = <9 M + igA^ the covariant derivative (in covariant gauges, ghost operators 
may also contribute in general, but not to physical matrix elements). 

These operators occur in the expansion of the product of electromagnetic currents at short 
distances, 

00 

J^\x) J^(0) = EE C nJ (* 2 , /x 2 , a s (^ 2 )) x Ml . . .x, n Oj^(0) , (119) 

n=0 / 

where they are distinguished by the behavior of their coefficient functions near the light-cone 

x 2 = 0, 

C nJ (x 2 , /i 2 , a s (fi 2 )) ~ (x 2 )- 2 hj(x 2 fi 2 , a s (/i 2 )) . (120) 

This singularity is identified by dimensional counting: the product of two currents has (mass) 
dimension 6, while the operators in (|117 ) and (|118|) have dimension 3 + (s — I) — 2 + s, with 
s the (maximum) spin (i.e., number of vector indices), while the corresponding factors of x^ 
contribute mass dimension —s. The power behavior of the coefficient function for any such 

5 Notice that d/dy~ = d + = <3_. 
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tensor operator of spin s and dimension D is thus (x 2 ) 3+( - D s ^ 2 . The short-distance expansion 
organized according to light-cone singularities is known as the light-cone expansion f42fl . 

The quantity D — s is called the "twist". All operators in eqs. ( p. 171 ) an d (|118|) have twist 
equal to two. Twist controls singularities on the light cone, and hence the high-g 2 behavior of 
the Fourier transforms of the products of currents, the DIS structure function. For DIS, then, 
the effect of Minkowski space is to elevate an infinite set of (twist-two) operators to leading 
behavior. 

Note that the light-cone x 2 = corresponds to the manifold x + = when the momentum 
p is in the plus direction. It is thus not the light cone along which the target particle moves, 
but rather the opposite-moving light cone, corresponding to a light-like scattered parton in the 
"brick-wall" frame. Of course, the scattered parton is not always in the opposite direction, but 
all the details of final states are absorbed into the hard scattering function H of fig. pj] and 
eq. (^). This is another consequence of factorization. From the point of view of calculating 
long-distance behavior in the parton distribution, the entire scattering process may be replaced 
by a pair of Wilson lines, as in eq. (|9l"l). There is a strong similarity here to the effective 



field theory picture often used to discuss the dynamics of heavy quarks f43fl , and indeed, the 
symmetries of heavy-quark effective theory are closely related to the universality properties of 
parton distributions, to which we now turn. 



4.6 Hard hadron-hadron scattering 

Once we introduce the concept of factorization, it is natural to apply it beyond inclusive DIS 
Jl|, |44j]. Of course, we must be careful to consider only inclusive hard-scattering processes, for 



which we may hope to find the necessary incoherence between the short-distance scattering and 
long-distance hadronic binding effects. 

In an important example, we consider processes in which a quark from hadron h annihi- 
lates an antiquark from hadron h', forming a virtual electroweak vector boson, which decays to 
a lepton pair. This reaction, 

h(p) + h'{p') -> l£'(Q») + X , (121) 
with its characteristic signal of a lepton pair (momentum Q M ), is known as the Drell-Yan process 



4"5fl . Its observation was one of the early successes of parton ideas, especially because it signals 
the presence of of quark pairs within ordinary hadrons. 

It is easy to write a factorization formula for such a cross section, by a straightforward 
generalization of the expressions for DIS structure functions, eq. (PU|). It is a convolution of 
two parton distributions, one for each hadron, with a hard-scattering function H. At lowest 
order, H is given by the Born cross section for quark pair annihilation to the relevant leptons. 
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At higher orders, gluon-quark scattering may also contribute, 
da hhl ^ Q2 (s,Q 2 ) _ ^ f 1 JtjH ± (t _ 2NtT ( Q 2 Q 2 



dQ 2 



J2 I' d£dg 4>i/h{L^)H i: 

n ~Z J 



i,j=f,f,G 



« s w ( 122 ) 



As in DIS, the hard-scattering function is a power series in a s (fi 2 ). H depends on the scheme 
chosen for the parton distributions. As an example, for -ff//, we have, to one loop in DIS scheme 

HI, 



H 
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, (Born) - 
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ln(l - z) 



+ 3 



1 



—6 — 4z — In 2) + 



4vr 2 



l) *(1 -*)]}) 



(123) 



where z = Q 2 /££'s. Given phenomenological parton distributions in some scheme, the factor- 
ization formula gives an absolute prediction for the Drell-Yan cross section, which has been 
successfully applied to a wide range of experiments. The corrections in H are not always small, 
however, and as we shall see, we sometimes need information about contributions at arbitrarily 
high power. 

Another application of parton model ideas, extended to perturbative QCD, involves single- 
particle inclusive cross sections, which count hadrons at fixed momenta, but are otherwise 
inclusive in the hadronic final state, 



h(p) + ti(p')^C(p c )+X. 



(124) 



If the hadron (C) is observed, for instance, at large transverse momentum, we know that a 
hard scattering has taken place, and may hope that incoherence and hence factorization is 



relevant [46], |47| . In this case, the parton model suggests that the hadron C arises from the 
"hadronization" , or fragmentation, of some parton k. The process of hadronization should, 
following our discussion of Section 1, occur over time scales that are independent of the hard- 
scattering scale, and of the fragmentation of other partons, scattered in other directions. Hadron 
C is thus expected to be produced in a universal fashion from parton k, and to inherit a 
fraction < z < 1 of that parton's momentum. The (incoherent) probability for this evolution 
is summarized in a "fragmentation function" dc/k{ z , fJ 2 ), which describes the distribution of 
hadrons in the fragments of a parton, and is analogous to the parton distribution <fii/h, but 
with the roles of hadron and parton reversed. In perturbation theory, d must be renormalized, 
and thus it depends on the factorization scale /1. The corresponding factorization formula for 
single-particle inclusive cross sections is 



uj c - 



da, 



hh'-+C(p c ) 



(P,P',Pc 



d 3 p, 



c 



E 



i,j,k=f,f,G 







1 dz 
d^'—H ilk 

z 1 



fi 2 pc-ip Pc-Cp' ( 2^ 



x 4>i/h(£,n 2 )4>j/h>(£,n 2 ) d c/k (z,fi 2 ) . 
45 



(125) 



The extra factor 1/z 2 allows H to be normalized to the corresponding Born cross section for 
i + j -> k + X. 

Like parton distributions, fragmentation functions are universal, within a given scheme to 
define them, and the same functions appear in hadron-hadron scattering, DIS and e + e~ anni- 
hilation. Also like parton distributions, their //-dependence may be analyzed, and summarized 
by evolution equations |39fl . We shall not go into these applications here, however. Rather, we 
shall close this section with a few comments on how generalizations of DIS factorization are 
established in perturbation theory. 



4.7 Jet-soft analysis 

The proof of factorization theorems fl], f|4], ^8] like those described above is highly nontrivial 



in perturbation theory, and, indeed, has not reached the sophistication of technical treatments 
of the operator product expansion in Euclidean space. We may, however, briefly discuss a few 
relevant physical issues. 

The essential complication in demonstrating the factorization of hadron-hadron cross 
sections is evident in the relevant pinch surfaces, shown in fig. The jets of the two incoming 



p 




Figure 20: General reduced diagram for pinch surface of a Drell-Yan cross section. The subdi- 
agram H includes possible final-state jets, and for lack of space the Drell-Yan pair is reversed 
in direction. 
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hadrons are connected by soft gluons, whose momenta vanish at the pinch surface. Roughly, this 
corresponds to the scattering of partons in each hadron from the Coulomb fields of partons in the 
other hadron before and/or after the hard scattering. Final-state interactions are present as well 
in DIS, and we expect them to cancel by unitarity arguments; initial-state interactions, however, 
are special to hadron-hadron scattering, and could, in principle, lead to a rearrangement of 
color, and even transverse momentum, of the partons in each hadron, due to the presence of 
the other hadron. Such a rearrangement could, in turn, affect the hard scattering, and break 
the universality necessary to identify the parton distributions of hadron-hadron collisions with 
those of DIS. This does not happen, however, for the following physical reasons. 

We recognize that the quanta of perturbation theory are gauge fields which include 
unphysical as well as physical degrees of freedom. Physical information is associated with field 
strengths, F^, and the latter behave very differently than the former under Lorentz boosts. In 
particular, the field strengths are strongly Lorentz contracted, while the unphysical polarization 
of A^ 4 , proportional to the momentum of each gluon, actually grows under Lorentz boosts. 
Thus, we expect the effects unphysical polarizations to contribute in individual diagrams in 
perturbation theory, but to cancel in a gauge invariant sum over diagrams. 

Such an analysis requires, in addition to the identification of pinch surfaces and power 
counting, an application of the Ward identities of the theory that decouple unphysical gluons 
from physical processes. Let us sketch how this can be done, taking as an example connections 



of soft lines to a final-state jet, as shown in fig. [H. 



(a) 



O 



H 






Figure 21: (a) Reduced diagram in which soft lines connect to a final-state jet. (b) Factorization 
of soft lines. 
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At a typical pinch surface, soft lines, q to the left of the cut, and q' to the right, are 
attached to the jet, which we may take in the plus direction, with large momentum component 
p + . We may define both soft momenta to flow into the jet, and then along jet lines away from 
the final state and toward the hard scattering functions H. For each jet line along which a 
soft momentum flows the momentum is then of the form i — q for lines to the left of the cut 
(amplitude) and t + q to the right (complex conjugate amplitude). Here £ is the momentum of 
a jet line in the absence of the extra soft momentum, and is hence naturally parameterized as 

E = (x P + ,r,e ± ), (126) 

with < x < 1. For I a jet momentum, we can assume that xp + is its largest component. We 
can now make two approximations that will enable us to factor the soft gluons from the jet, 
replacing their couplings by an effective eikonal line. The first is that the coupling of the gluons' 
propagators G u \(q) to the jet is always through the combination p u G u \(q), that is, that the 
soft gluons couple only to the large component of the current. The second is that propagator 
denominators (I + q) 2 , depend only upon the minus component of the soft gluon's momenta, 
i.e. the components in the direction opposite to the jet direction, 

(£ ± q) 2 ~ ±2xp + q- + f . (127) 

We shall return to this condition in a moment. 

Once these approximations are made, the jet depends only on one component of each 
soft gluons' momentum, and on the same component of its polarization. Thus the coupling 
of soft gluons to a jet moving in the plus direction is equivalent to the coupling of a set of 
gluons with only minus momenta and only minus polarizations. Such gluons are unphysical, 
and are equivalent to a phase rotation on the external lines of the jet, which attach to the 
hard scatterings. In an abelian theory, for instance, an arbitrary jet with one vector (photon) 
attached to each side of the cut, is equivalent to a jet with no photons, multiplied by two 
linearized propagators ("eikonal lines"), one for each of the photons, 

r/(p, q, q') ~ J+ + (p, q, q') = g 2 - f , / J (p) , (128) 
where Jo(p) is the jet with no soft photon connections and where 

r = (0,q~,0). (129) 



The factors on the right of ( 128 ) are of the form of the propagator and vertex in eq. ( p2|) with 



n = p. The eikonal lines here are illustrated in fig. plfb. 

In a nonabelian theory, multiple gluons still factor, and their color interactions are sum- 
marized by connections to eikonal Wilson lines, which have the same perturbation theory rules 
for lines and vertices as those in eq. (|9~i|) for DIS parton distributions. Schematically, we write 
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where the E's are generated from the operators 



U(A) = exp 



-w 



dX p ■ A{Xp) 



(131) 



[TJ, due to the unitarity of the Wilson 



The soft divergences factorized in this manner cancel 
lines, 

U[A]* U[A] = I . (132) 

To realize this cancellation, it is necessary to sum over all final states that differ by soft gluon 
emission. This, of course, is exactly what we do in inclusive hard-scattering cross sections. 
Given this cancellation, the remaining jets are independent of each other, and may be factored 
from the hard scattering as in DIS. 

The argument that we have given above is, of course, quite rough. In particular, the 
approximation ( 127 ), is highly nontrivial It is clearly not true for every soft momentum g M , 
since it fails whenever q~ vanishes compared to the other components of g M . The approximation 
will hold, however, if the q~ integral is not trapped at or very near zero. Arguments to this 
effect are relatively easy to give for may cross sections in e + e~ annihilation, where we verify 
that all the poles in q~ from lines within the jet, (£ ± q) 2 =F ie, are in the upper half-plane 
for soft gluon momenta routed as above. The situation is much more difficult when the jet 
originates from an incoming hadron, as in the Drell-Yan process. Here, pinches may occur on 
a diagram-by-diagram basis, as in fig. where the lines xp + q and (1 — x)p — q have q~ poles 
in opposite half-planes. These pinches cancel, however, in sufficiently inclusive cross sections, 



(l-x)p-q — ► 








xp+q 




q 





Figure 22: Reduced diagram for which q in pinched close to zero. 

for the physical reasons that we have outlined above, and soft gluons that are not in the jet 
direction decouple from incoming as well as outgoing jets |1|, |48| . 



We now ask when a soft gluon is close enough to a jet's direction to be part of that jet. 
Let us consider the case when the two momenta in eq. (|127|) are both on-shell, i.e., q 2 = £ 2 = 0. 
To factor q from the p-jet we require, in particular, that 



xp + q > \£± ■ q±_\ 



(133) 



This condition has a natural geometrical interpretation, which itself has important physical 
consequences. For a well-collimated jet, 0£ = \£±\/(xp + ) is the angle of £^ to the overall jet 
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direction. Defining 9 q similarly, and parameterizing q + = yp + , we have 

\t±\ = e exp + , \q±\ = e q y P + , <r = ^|+. 

Using these relations it is easy to show that the soft gluon q will factorize from the jet, i.e. eq. 
( |127| ) will hold, unless 

9 q <9i, (135) 
that is, unless q is emitted closer to the jet axis than £, the particle that emitted it. This 



"angular ordering" |0J, limits the region of phase space into which soft gluons may be emitted 
late in jet evolution, and hence suppresses the multiplicity of very soft gluons, since gluons that 
factor from the jet do not get enhancements from pinch surfaces where jet lines are near the 
mass shell. This is one of the many consequences of the coherence of gluon emission in QCD 



5 Two-scale problems I: Sudakov resummation 

The factorization theorems of the previous section go far toward connecting perturbative QCD 
to experiment. Their application, however, is limited somewhat by our assumption that there 
is only a single large scale in the problem, for instance, Q 2 in DIS. In this case, it is natural 
to choose the factorization scale to be of order of Q, and to use the evolution of the parton 
distributions discussed above. The coefficient function is then of the general form 

C{xl£) = Y,d n {xl€)(* n a , (136) 

n 

where d n (z) is a distribution in z. Then, if a s (Q 2 ) is small, we may hope that the effect of 
higher-order terms in C is small as well. 

Unfortunately, this is not always the case even in DIS; the d n (z) will generally include 
factors like a™ ln n (z)/z, or even a™ ln 2n_1 (l — z)/(l — z). When x is very small, or very 
close to unity, these logarithms may produce large corrections at every order. Recalling that 
(p + q) 2 /Q 2 = (1 — x)/x, we see that in both cases the arguments of the logarithms are ratios 
of kinematic invariants. There are many other such examples, all with two hard scales, Q\ and 
Q 2 ,, which satisfy 

Ql > Ql > A 2 . (137) 

The first inequality ensures the presence of a large ratio in the hard scattering function. The 
second ensures that parton densities evolve perturbatively, and that the running coupling in 
the hard scattering function, a s (Q?,) is small, even if the combination a s {Q^) \n{Q\/Q^) is not. 
In this case, we may try to "resum" the series in \n{Q\/Q^) to all orders in perturbation theory. 

We shall treat three representative cases: T — > 1, with T the thrust, eq. (|80|), in e + e~ 
annihilation, Qt <C Q for Drell-Yan production of pairs of mass Q, and x — > in deeply 
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inelastic scattering. The first and second illustrate the resummation of so-called "Sudakov" 
double logarithms, and the second and third illustrate resummations based on an extension of 
factorization to fixed transverse momenta. 



5.1 Sudakov double logarithms 

Let us return for a moment to our dimensionally-regulated expression for one-loop corrections, 
eq. (fUD, to the electromagnetic form factor, 



T/^ 2 , e) = -«e// u (pi)7^(p 2 ) p(q 2 , e) . (138) 

At one loop, its momentum-dependence is contained in an overall factor (— q 2 )~ e , which we may 
expand to order e 2 to get an expansion in lng 2 , 

p( q \e) = l-^-C F ln 2 (q 2 /^) + .... (139) 

The double logarithms in momentum transfer are another reflection of the overlap of collinear 
and soft singularities in the vertex function, and are generally referred to as Sudakov dou- 



ble logarithms [f| [50], [51]. Although the electromagnetic vertex is not infrared safe, it is an 
important element in any process whose kinematics approaches those of elastic scattering. 

Our first example |52 is the e + e~ annihilation cross section near unit thrust, T — 1. Here 



the double logarithms appear not in Q 2 directly, but in (1 — T). Consulting the definition of 
thrust, eq. (|8~0"D, we see that at T = 1 the final state consists of two back-to-back massless jets. 
A little kinematics also shows that when both jet masses pf are small compared to Q 2 they are 
related to the thrust by 

1-T = ^i±^. (140) 

Q 

At order a™ the leading logarithm in 1 — T is given by an exponential of Sudakov loga- 
rithms, 

1 da -2CV^ ln(1 - T) exp i-CA ln 2 (l -T)\ . (141) 



Otot dT tt 1 - T { 7i 

As T approaches one, ^ vanishes. This is a quantum-mechanical reflection of the classical 
radiation field that must accompany any process in which a charged particle is accelerated 
(quark pair creation being an extreme example). Quantum field theory assembles this classical 
field out of many soft and collinear gluons (the correspondence principle). Cross sections in 
which gluon emission is forbidden in part of phase space are correspondingly suppressed. 



5.2 Factorization for T — ► 1 



Our goal in the following is to rederive eq. ([141 ), and to extend it to include nonleading 
logarithms and the effects of the running coupling. Our arguments apply to a large class of 
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cross sections. In fact, the resummation of Sudakov logarithms follows from the factorization 
properties of the cross section in the regions of momentum space that give rise to the logarithms. 
This factorization has already been illustrated by the physical picture of annihlation given in 



fig. |T^, which shows a two-jet configuration. The infrared divergences associated with such 
configurations cancel, according to the arguments of Section 3.3, but for T ~ 1, the cancellation 
between diagrams with virtual and real gluons is constrained by the requirements that the real 
gluons be either very soft or emitted very close to the quark or antiquark directions. This leaves 
large finite remainders in the cancellation of divergences. These are the logarithms of 1 — T. 

Let us recall the power counting arguments of section 2.3. We saw there that in an 
axial gauge, £ • A = 0, all collinear divergences may be absorbed into the jets. As a result, 
in such gauges, double logarithms arise from collinear gluons in jets which become soft, while 



staying collinear. The form of fig. [IJ3, combined with the factorization of soft gluons from jets, 
described in Section 4.7 above, suggests that the thrust cross section, like Drell-Yan, factorizes 
into functions that describe the two jets, the hard, and the soft subdiagrams. At double 
logarithmic accuracy in axial gauge, therefore, we may neglect soft gluons that connect the 
jets to each other. This approximation will simplify the arguments below, without changing 
the character of a more general treatment, which gives an essentially equivalent result, but is 
accurate to all logarithmic order. At fixed values of jet masses p], we shall therefore begin with 
the factorized expression 

dp 2^2 = J l(Pl>M)^2(p2,^,0#(PbP2,/i,0 , (142) 

where the J's represent the jets, and H the hard-scattering factor. We shall suppress particle 
labels, but it is relatively easy to show from power counting that to leading power in 1 — T the 
jets are connected to the hard scattering by a quark and antiquark only. Eq. ( |142|) is enough 



to derive the exponentiated double logarithms of eq. (|141|) above, with corrections due to the 
running of the coupling. 

In axial gauge the jets depend, not only on their invariant momenta, but also on their 
energies through the products Pi The jets and H are thus not individually Lorentz invariant. 
For a general gauge vector the precise arguments are Pi-C./V^ 2 because the gluon propagator, 
and hence J is invariant under simple rescalings of (In the following, we set £ 2 = 1.) 

In view of eqs. ( |140| ) and ( |142|) , the cross section at fixed 1 — T ~ is of a convolution 

form 

~ - J o Q2 dpldplS^-T-^-^H(p 1 ^/^p 2 ^/^a s ^)) 

x Ji (pi//i 2 ,pi -£//i,a s (/i 2 )) J 2 (pI/^ 2 ,P2 -^//i,a s (/i 2 )) . (143) 

Dividing by the Born cross section o"o gives H = 1 at lowest order. Because we are interested 
in the limit pf/Q 2 — > 0, the Pi may be expanded about back-to-back lightlike momenta. We 
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choose axes so that Pi is in the plus direction and p 2 is in the minus direction. Then pf and 
p 2 are both nearly equal to Q/\/2, while 

(144) 

For £ in an arbitrary direction, we may take p\ • £ ~ and £>2 ' £ ~ P2 £ + to leading power 
in Q. Then to leading power in 1 — T, the pf integrals are independent of pi ■ £. 

We want to identify singular ln m (l — T)/(l — T) behavior for T — >• 1, and for this purpose 
moments with respect to T are particularly useful, 

~a{n) = - f dTT n ^, (145) 

since the moments of any function that is finite at T = 1 falls off as 1/n for n —>■ 00. In 
particular, logarithms of 1 — T are transformed into logarithms of n by 

1 t n — 1 C_i ") m 

dT — ln m (l — T) = i ^- ln m+1 n + . . . . (146) 

l— 1 m+1 

Keeping only terms that are finite or grow as n —> 00, and using the relation T n ~ e~ n( - 1 ~ T \ 
which holds in this approximation, we find that the convolution in (|143|) factors into a simple 
product under moments, 

a(n) = Ji (Q 2 /n{i 2 ,p l ■ ^//i, a s (/i 2 )) J 2 (Q 2 /n^ 2 ,p 2 ■ ^//i,a s (/i 2 ; 
*H fa ■ £/n,p 2 ■ ^//i,a s (/i 2 ; 

up to corrections that vanish as 1/n. Note that n now appears only in the combination Q 2 /n/i 2 . 
By analogy to the derivation of evolution for DIS structure functions, we shall use this factorized 
expression, coupled with renormalization group arguments, to derive a resummed cross section. 
The new feature in our Sudakov factorization is the dependence on the axial gauge vector £. 
Although each of the factors that makes up a is ^-dependent, the physical quantity must be 
gauge- invariant. This invariance will drive the resummation pT| . 



5.3 Resummation for T — > 1 

We start with the renormalization group behavior of of J, which is simple, since it has only 
two external (quark or antiquark) lines, 

lnJ=-2 7? , (147) 

with 7q the quark anomalous dimension. Since the cross section is independent of the renor- 
malization scale /i, H must behave in a corresponding fashion 

\nH = 4 lq . (148) 



— + P— 

^ dp, dg 



— + 3— 
^ dp dg 
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We now observe that very similar reasoning may be applied to the gauge-fixing vector £. 

Let us change (boost) £ + and £~ in a manner that leaves £ 2 = 1. Using the discussion after 
eq. ( |142| ), the independence of cr(n) from £ M may be expressed as da/dhi^ = —da/dln^ + . 
The chain rule then gives 

d\nJ 2 (Q 2 /n/j 2 } p 2 //j,) | 8\nH(pt/n,p2 /» _ d In J^Q 2 jn\x 2 ,pf / ' jj) d In H(pf/n, p£jjj) 
dlnp 2 d\np 2 dlnpi dlnpf 

(149) 

where we have made components explicit and have suppressed a s . This relation is surprisingly 
powerful, because Ji, J 2 and H depend on different sets of arguments. J 2 , for instance, depends 
on both p 2 and Q 2 /nfi 2 . Its derivative with respect to p 2 may depend upon either of these 
arguments, but must be cancelled in (|149|) by the derivatives of H and Ji, which contribute 
additively. Thus, its dependence on p 2 and Q 2 /nfi 2 can only be additive after the derivative, 

9 \nJ 2 (^-^,a s ^ 2 )) =K(^a s ^ 2 ))+G(^a s ^ 2 )) . (150) 



The function .K" cancels a corresponding term from Ji while G cancels the contribution from 
H, whose derivatives must satisfy, 

dlnp'" dlnpi ^(/i) ^ ( /i ) ' (151) 

This separation of short-distance and long-distance dependence in jets is characteristic of Su- 
dakov factorization. As the gauge changes, the jets exchange contributions with each other (via 
K) and with the hard part (via G). Here we find a strong analogy to the "matching conditions" 
of effective field theory. 

We now have at our disposal two evolution equations, the first relying on invariance 
under renormalization group rescalings, the other on gauge invariance, but both based on 
factorization. Combining the two, we shall find enough information to determine all logarithmic 
n-dependence. 

By eq. (p.47| ), (dlnJi/d/j,) is independent of momenta, so, for instance 

—^-lnJ^O. (152) 
a/j dpi 

Applying this result to ( |150|) , we conclude that the combination K+G is itself a renormalization 
group invariant 

IX— (K + G) = 0, (153) 



which implies that yet another anomalous dimension relates K and G [5(| [51], p3 |, 

/I— K = - 7fe (cO = -/i—G . (154) 
d[i d[i 
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Now we can relate the moment-dependence of K and G through 

Q 



K ( -^,a s (/i 2 ) 



G I — ,a s (/i 2 ) 
V n 



K(l,a s (Q 2 /n))+G(l,a s (Q 
1 fQ 2 , /2 , 



(155) 



1 fQ 2 dfjf- 

IQ 2 /n // 

in which all logarithms of n are generated either through the running coupling and/or the 
explicit fi' integral. There are only two steps left, to solve for the n-dependence of the J's and 
to combine everything together in the cross section. 

Combining eqs. (|150|) and ( |155| ), we derive the full evolution of J2 in terms of P2 , an 
exactly similar equation holding for J\ in terms of pf, 

1 rQ 2 d\ 2 







■ In J 2 = K(l, a s (Q 2 )) + G(l, a s (Q 2 )) - - / — Tj(a s (X 2 )) . 

2 JQ 2 /n X Z 



(156) 
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Here Tj combines 7^ with a term that allow us to have the same running coupling in K and G, 



d 

Tj(a s ) = j K (a s ) + p(g)—K(l,a s ) 
— )2C F +[ — ^ 

7T / V 7T 



(157) 



In the second line, we have given the two-loop expression for Tj, where as usual n/ is the 
number of quark flavors. 

It is now a simple applications of the chain rule to derive a differential equation for the 
n-dependence of both jets, 







-el 



In J (^1, % a , ( A = 1 rs („.(„»)) - 1 f (q,(A 2 ; 



9 In n 2 dg 
where we have set pi ■ £ = Q and where 

T'j (a s (u 2 )) = Gj (l, a s (Q 2 )) + fG (l, a s (Q 2 )) - 2 7? (e* s (u 
The solution of eq. (|158| ) relates J at large n to J at n = 1 



158) 



(159) 



In J — -, — , a s (n ) 
\njjL z fi 



In J — , — ,a s (n n 
1 rQ 2 d\ 2 l u 



2 jQ2/ n X 2 



ln^r, {a s (X 2 ))-Tj' (a s (X 2 ] 



. (160) 



If we now set /1 = Q, all logarithms of n are generated by the integrals of the two anomalous 
dimensions Y and V and the expansion of a s (Q 2 /n), and, as promised, exponentiate in the 
moments of J. 



6 One way to verify this result is to observe that (3(g)dT (a s (X 2 )) / dg = XdT(a s (X 2 ))/d\. 
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The resummed expression ( |160| ) for the jets organizes all logarithms of n in the moments 
of the cross section, since the hard function H has no Inn-dependence. The inverse transform 
of J(n) then gives the singular 1 — T-dependence. To be explicit, the inverse Mellin transform 
of J in ( |160| ) is given by [[52], [53[ 

- e E(i-x,a s (Q*)) n sin ( vri p 1 )r(i + p x ) + . . 



J 1(1 - x)Q-, — ,a s (u 2 ) 

V hp 



\l*=Q- 



1 — X 



(161) 



where 1 — x = p 2 /Q 2 , and where the exponent E is given by the right-hand side of (|160|) with 

fi = Q, 



E 



'(l-x)Q 2 A 2 

while Pi is related to the exponent by 

Pi = 



— [ ln(^)r J (« s (A 2 )) - r'j(a s (A 2 )) ] + In J (l, 1, a s ((l - x)Q 2 



dE(l-x } a s (Q 2 )) 



(162) 



(163) 



dln(l — x) 

Terms omitted in eq. ( |161| ) are suppressed by powers of p 2 /Q 2 . It is a straightforward matter 
to verify that the leading logarithms in 1 — T in eq. ( |141| ) are indeed generated by this form. 
Here, however, we see the essential role of the running coupling for nonleading logarithms. In 
particular, because of asymptotic freedom, the exponent receives relatively larger contributions 
from long distances, and relatively smaller contributions from short distances than in the case 
of a fixed coupling. We shall return to the consequences of this observation in the final section 
below. 



5.4 kr- factorization for the Drell-Yan cross section 

As another example of the variety of interesting cross sections to which a variant of ( |142j ) 
applies, consider Drell-Yan cross sections at measured pair mass squared Q 2 and transverse 
momentum, Qt- (The discussion below follows the extraordinary analysis of Collins and Soper 
for transverse momentum distributions in e + e~ annihilation [pl[|.) Here again a factorization 
holds in covolution form, but now the convolution is in terms of the transverse momenta of 
gluons emitted from jet functions associated with the incoming hadrons, along with "central" 
soft gluons from the soft subdiagram of fig. For the Sudakov resummation of logarithms 
of Qt in the Drell-Yan cross section, a convolution in transverse momentum will play the role 
of the convolution in jet mass for the thrust distribution. Otherwise, the reasoning is quite 
similar. 



Explicitly, the convolution is |55 



ddhh' \- { AC A? 1 f 

dQ 2 d 2 Q T ~ 2fJ ( * < * J (2tt)6 



xP f/h & kr) P f/h ,{i', k' T ) H ff (Q 2 ) S(k T , s )5 2 (Q T -k T -k' T - k s , T ) , (164) 
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which factorizes under a Fourier transform, 

- 'Ljj^Plll>{t^,QIAP} /h .{(:,b l i,Ql l i)S(b l i), (165) 

where <To = Aira 2 /9Q 2 s, the Born cross section, summarizes the hard part H to leading order 
in a s . The jet functions are defined as matrix elements of quark fields separated by a spacelike 
vector (0 + y _ , b), 

Ps/ h (^b^p + /ri = I ^-e-^- < h(p) | g / (0 + ,y-,b) 7 + g / (0) | h(p) > , (166) 



47T 

where we suppress an average over spin. These matrix elements are gauge-dependent, but in 
any axial gauge they absorb all double logarithms of b (or Qt in momentum space). At b = 0, 
they are normalized to the quark distributions eq. (|88|) . The gauge vector plays the same role 
as in the case of thrust, and the jet matrix elements obey a noncovariant evolution equation 
that separates their dependence on hard and soft scales, 

9 , » rs U t 2x\ , „ (Q 



lnP = Kp(bfi,a s (fi 2 ))+Gpl^,a s ( f i 2 ) \ , (167) 



dlnQ 

where the combination Kp + Gp is renormalization scale invariant, 

H^-{Kp + G P ) = Q. (168) 
djj, 

Following essentially the same reasoning as for the n dependence of the jets in the thrust cross 
section, the logarithmic dependence of the jets exponentiates in the transform (b) space, 

1 rQ 2 d\ 2 



exp 



i/b 2 A 2 




<j> f/h (^l/b 2 )(l + a s (l/b 2 )), 

(169) 

with Tj as in eq. ( |157| ) and B a power series in a s . This approximation holds in the range of 
b for which Q » 1/6 >> A, and organizes all perturbative logarithms of Qt/Q in the cross 
section. 

Sudakov resummation may be relevant to any cross section with an underlying hard-soft- 
jet factorization. The exponentiation of logarithms requires a convolution in phase space, like 
( 164|) or ( |143| ). Applications include threshold corrections (where the relevant variable is 1 — z = 
1 — Q 2 s) for the inclusive Drell-Yan cross section |54|, |56| (|122|) as well as other, purely QCD 
cross sections such as top or jet production |^7|. Another important case involves transverse 
momentum distributions in e + e~ annihilation (where the first really complete analysis of such 



a process was carried out [51]). Yet another example is semileptonic B meson decay at the 



endpoint of the lepton energy spectrum [58], where the lepton recoils against a jet of hadrons. 
Undoubtedly, there are others as well. 

We now turn to another classic resummation of large logarithms, organized by the BFKL 
equation. 
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6 Two-Scale Problems II: Small x and the BFKL equa- 
tion 



The small- a; limit of deeply inelastic scattering is one of many cross sections that show a set 
of enhancements organized by the BFKL (Balitskii-Fadin-Kuraev-Lipatov) equation || [59]. 



Here again a transverse momentum factorization may be used as a starting point, although 
in a rather different kinematic region from the Drell-Yan cross section just discussed. In DIS, 
these enhancements appear as logarithms in x at fixed Q 2 . For Q 2 not very large, the quantity 
a s (Q 2 ) ln(l/x) can be large, and we may be tempted to resum corrections of this sort to all 
orders. We begin with a brief review of the origin of logs of x in the evolution formalism 
developed in Section 4 above. 

6.1 x -> for DGLAP evolution 

Referring to eq. ( |101|) above, the kernel PqJ g which describes gluon-to-gluon evolution is sin- 
gular as x — > 0, so that in the notation of ( |115| ) (and using Ca = N), 

, . 2Na s 1 

Pgg(x) = - + .... 170 

71 X 

This behavior produces a pole at n — 1 in the corresponding diagonal element of the singlet 
anomalous dimension matrix, eq. ( |114| ), 

1gg(ji) = + . . . . (171) 

7r n — 1 

We now recall the solution to the renormalization group equation for the moment of a 
DIS structure function F(n,Q 2 ), which follows from eq. ( |112| ). 

F(n, Q 2 ) = F(n, Q 2 ) e ^ )/b * )lnt , (172) 

where ln (a s ) = ^ (a,/ir) + . . ., and t = HQ 2 /A 2 )/ HQ 2 /K 2 ). 

For x — > 0, the (inverse) transform from F(n, Q 2 ) to F(x, Q 2 ), 

F(x, Q 2 ) = — e -^ X +4N/(b 2 (n-l)) Int ^ (1?3) 

J-ioo 2rri 



has a sharp saddle point at (n — 1) = J4N\nt/ (62ln(l/x)), which gives the x behavior 



F(x,Q 2 )~ eVW^MWw. (174) 

This striking result shows a rapid increase as x — > 0. It relies, however, on DGLAP evolution, 
which assumes that In Q 2 is relatively large. It is natural to ask what happens if x is so small 
that, for instance, ln(l/x) ^> ln(Q 2 /A 2 ), and to treat the resummation of logarithms of x 
self-consistently. 
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6.2 /^-factorization for DIS; the BFKL equation 



The standard DIS factorization, eq. (0), assumes an ordering in transverse momenta, which 
allows us to decouple the transverse momenta of the partons from the hard scattering. Now, 
however, we do not wish to treat Q 2 as arbitrarily large, so we generalize (|90|) to a convolution 



(175) 



in both £ and kr- The fcr-factorized form of a DIS structure function F(x, Q 2 ) is |]60| , |B~T 

F(x, Q 2 ) = J d 2 k T (~ Q, k T ) ^(£, hr) , 



with hr) a generalized parton distribution at measured x and k T , and C(x/£,Q,kr) the 
corresponding coefficient function. Since leading logarithms of x are generated by purely gluonic 
evolution, we shall restrict ourselves to gluon distributions, neglect mixing, and suppress parton 
indices. This factorization is illustrated in fig. EBI. It is appropriate for the limit Q 2 fixed, x — > 0. 





F 



Figure 23: /^-factorization in DIS. 



The leading logarithms of x in ( |1 75|) are generated by a large set of diagrams. The 
diagrams that show the relevant mechanism most clearly are the ladders, illustrated by fig. 



To understand how logarithms of x are produced, we parameterize the momenta flowing on the 
sides (vertical lines in the figure) in components parallel to the incoming momenta p and q and 
transverse components, 

ki = dip + (3iq + k iT . (176) 

Here q is approximated by a lightlike vector, because as x — > the ratio of the photon invariant 
mass to the center of mass energy vanishes. Logarithms of x result from configurations in which 
the "light-cone" fractions q.{ and j3j are strongly ordered, but the transverse momenta are all 
of the same order, 



«i >a 2 3> • ■ ■ 3> 

Pi < fa < ' ■ ■ < Pn-1 
kiT ~ kj T ■ 



(177) 
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n - 1 




Figure 24: n-gluon cut ladder. In DIS, the photon would attach to a quark loop at the top. 

It is the lack of ordering in transverse momentum that distinguishes BFKL evolution. Useful 
approximations that follow from ( |177 ) are 



£A 

1=1 

n-1 
i=j 



0. 



J ! 



~ a 



for the lightcone fractions and 



k ~ -ki 



LT 5 



178) 



(179) 



for the invariant masses of lines on the sides of the ladder. The emitted "rungs" of the ladder 
carry momenta 



ki—i ki 



;i80) 



with k = p and k n = —q. 



A diagram like fig. with n rungs generates n — 2 logarithms of x, which come from 
ordered logarithmic integrals over the a«, or equivalently rapidities yi of the ith rung, 



Vi 



1 ln Oij-l 

2 n A 



;i8i) 



It is relatively easy to identify the leading behavior of the ladder diagrams. There are 
n mass-shell delta functions from the cut rungs, S(£^). The top-most mass-shell condition 
fixes a n _i ~ x, while the remaining 5(£f) fix n — 1 fractions i = 1 . . . n — 1. The leading 
numerator factor comes from the terms that contract the incoming momenta of the top and 
bottom vertices to give (p-q) 2 , while the vertices to which ith rung connects produce a factor of 
order \(k Tti _ 1 + k Tji ) -ej(£j)| 2 , with e^) the polarization of the ith emitted gluon. The resulting 
term is logarithmic in the remaining, ordered Qj integrals, which are each of the form, 



dcti/cti 



dyi 



(182) 



2/i-l 
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with a minimum value a n _2,min = x. Transverse integrals also give logarithmic power counting 
on a graph-by-graph basis. The resulting logarithms, however, cancel in the sum over diagrams. 
This cancellation will be evident in the equation we derive below. 

Consider what happens when we insert yet another gluon into an arbitrary diagram, as 
illustrated in fig. |25|. The diagram has been factorized because all of its gluons are already 
ordered as in ( 177 ). To get an additional logarithm from this gluon, it should fit into the 



-C 





Figure 25: Adding an extra strongly ordered gluon. 



rapidity gap between C and and should have a transverse momentum comparable to those 
of the other gluons. To get the logarithm, however, the "new" gluon, of momentum £, need 



not itself be inserted as a ladder. The ladder insertion, the first in fig. 2S[ works, but so do 
many insertions that connect the vertical gluon and either the top or bottom, or even some 
that connect the top and bottom directly. Similarly, there are additional possibilities for the 
insertion of a virtual gluon, also illustrated in fig. 

To determine when each of these diagramatic insertions can give logarithms we use strong 
ordering Ql77| ). Because of strong ordering, whenever the momentum £ flows along a line of 
momentum £' + £, say, in C, (£' + £) 2 ~ 2£'~£ + = 2aij3i>s, and analogously for attachments to 
J 7 , but with the roles of the fractions a and (3 reversed. 

Luckily, however, it is not really necessary to worry about each diagram individually. 
Instead, we appeal to the jet-soft analysis of Section 4.7 above. Strong ordering implies that C 
is sensitive to £ + only, and T is sensitive to £~ only. At the same time, it is clear that the + 
component of the polarization of £ is also dominant in its coupling to C, and the minus to T. 
We are thus in the situation of eq. ( |128| ), fig. [21], and as in that case, the sum of all attachments 
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Figure 26: Eikonal factorization of additional strongly ordered gluon. 



of the gluon I to either of the subdiagrams C and T factors from it as in fig. |26], with an eikonal 
factor 



-gQ 



abd ) 



Q ■ 



183) 



for C, and similarly for its connection to T. 

The result of attaching our extra gluon is then summarized by a kernel K(k', k), which 



we illustrate with a cut gluon in fig. |27|. The new gluon attaches at either end to a vertex 
represented by a circle. The eikonal factors themselves (|183|) are as usual represented by 
double lines. 

Because the integrals are ordered, the separation between T and C in eq. ( |175| ) is at 
a definite value of a, or equivalently of plus momentum £ + = ap + . We may thus identify the 



combination of J-(k) and K{k! , k) in fig. 27 with J-(k'), now evaluated at external momentum 
k! = k — i. This gives an integral equation for the /c-r-dependent jet function T . Since we know 
there is a logarithmic integral in £ + , we write 

Here a plays the role of a factorization scale, separating C and T. Thinking back to DGLAP 
evolution for ordered transverse momenta, we derive a new equation by taking a derivative with 
respect to a, 



J d 2 k T K{k' T ,k T ) T{i+,k T ) 



;i84) 



ol^-F(ol, k' T ) 



J d 2 k T K{k! T ,k T )T{a,k T ) 
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k' k , 



k' 



Figure 27: Graphical representation of the BFKL kernel for a cut gluon. 



d 2 ki 



7T 2 J (k T -k' T ) 2 



k' 2 

F(a, k T ) - ■^jF(a, k' T ) 



where we have defined 



k 



t 



rv r p 



185) 



(186) 



This is the BFKL equation, as it appears in DIS, and K(k' T , kx) is called the BFKL kernel. In 
the second line the kernel has been evaluated explicitly for the diagrams of fig. [27]. Details of 
the evaluation are given in Appendix C. This is only one of many forms in which the kernel is 
expressed, but it is one of the simplest. Clearly, the first term comes from real-gluon diagrams, 
for which kx is not identically equal to kx, while the second term comes from virtual gluons. 



6.3 Solution of the BFKL equation 

Because the BFKL equation ( |185|) is a convolution in kj<, and first order in the derivative with 
respect to In a, it is natural to express its solution as an expansion in functions that are powers 
in a and kx, 
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where Cu V and fi 2 are constants. We want to solve for ou, which determines the power-law 
a-dependence, and the associated power v for transverse momenta, v is to be considered as a 
complex number, so that iv is not necessarily imaginary. The extra —1/2 in the power of k\ 
is a matter of convenience. Together, iv — 1/2 serve as an anomalous dimension appropriate 
to evolution with a fixed coupling (see eq. (|112|) ). From the discussion of DGLAP evolution 
in Sec. 4.4 above, we recall that logarithms, and hence evolution, in parton distributions arise 
from transverse momentum integrals. The argument of the running coupling is thus naturally 
chosen as fcf>, which is effectively fixed in our case, because of the strong-ordered kinematics of 
eq. f77|). 



Substituting the ansatz ( |187| ) into the BFKL equation ( |185| ), we readily find an implicit 
expression that relates u and v (see Appendix C), 




i 



7T 




d z k T { ^ 



2 1 I I ,,2 J (t. U \2 \ ,,2 / QZ.2 (U Ul \2 



2k\{k T - k' q 



TJ 



The infrared divergences of each term on the right-hand side of this expression cancel, as may 
be seen by carrying out the integrals via dimensional regularization, and the explicit relation 
between ui and v is 

u( v ) = J^L (r^ (i v + I) _ ^(i)) , (189) 
where the special function ip is the logarithmic derivative of the gamma function, 

tjj(x) = — InT(x) , 
dx 

^(1) = -IE- (190) 
Another common form for the relation ( |189| ) is found by defining 

7 = ^ + 1, (191) 

in terms of which 

W ( 7 ) = (^(7) + ^(1 _ 7 ) _ 2^(1)) . (192) 

7T 

A general solution to the BFKL equation for DIS may be written as a superposition of 
power-law solutions, 

The small- a; limit is dominated by a saddle-point of u as a function of v, in much the same 
way as the small-x behavior from DGLAP evolution in eq. (|173|) above. In this case there is a 
saddle point at v = 0, as may be seen from the series found by expanding the ip functions, 

u ( u ) = (2 In 2 + f](-l) fc (2 2fe+1 - l)C(2ife + l)v 2k ^j , (194) 



F(x,k T ) = J dux-^ u U^\ . (193) 
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with Q{2k + 1) the zeta function ((x) = J2n(^/ nX )- This leads to the (famous) asymptotic 
behavior || 

F(x,q T ) ~ x^ 1 " 2 ^^)- 1 / 2 , (195) 

in which we see a power-law enhancement as x vanishes, much stronger than in DGLAP evo- 
lution, eq. QT7§. 

The BFKL formalism was originally developed to describe hadron-hadron scattering 



in QCD, both the total cross section and the closely-related Regge limit [62], t fixed, s — > oo. 



The t — Regge limit is related to the x — > limit in DIS by the following "translation" , 

(p + q) 2 = —Q 2 - s 

x 

F(x) -> a tot ~^°. (196) 

Thus, at t = 0, o"tot grows as a power of s in this approximation. In fact, the total cross 
section for hadron-hadron scattering does increase at high energy, although an uninterrupted 
power-law rise would violate unitarity, as embodied in the Froissart bound [JB^j . The BFKL 
behavior that we have just derived therefore cannot be the final answer. 

Although the dominant BFKL power (|195|) is infrared finite, in most applications, such as 
the total cross section, there is no natural perturbative scale at which to evaluate a s . The scale 
of kr introduced above is, after all, quite arbitrary. There are, however, two-jet correlations |)3 



for which BFKL resummation is naturally infrared safe, with a coupling fixed at a perturbative 
scale of the order of the transverse momenta of the jets. The dominance of resummed leading 
logarithms in such cross sections may require very high energy ||64|| . Another infrared safe 
application of the BFKL formalism is to hypothetical heavy-quark onium-onium scattering, 



in which the inverse size of the onium wave function serves as an infrared cutoff |65| . This 
model is serving as a valuable laboratory for the study of forward scattering in a self-consistent 
perturbative context. 

At least two sets of corrections can lead to a softening of the BFKL or "bare" pomeron, 
nonleading logarithms and nonleading powers (higher twist). "Nonleading logarithms" refers 
to higher powers of a s at a fixed power of Ins. One approach is to derive nonleading terms 



in the expansion of the kernel K(k,k') |]67|]. Another is to compute exchanges of multiple 



ladders f66j. Indeed, a multiladder generalization of the BFKL equation |6£| may be used as 



an inspiration for a picture of forward scattering in terms of two-dimensional field theories. In 
these investigations, QCD comes into contact with conformal field theory and the theory of 
exactly soluble models. 



In addition, ladders may interact with each other, a process that produces "shadowing' 
the softening of parton distributions due to the spatial overlap of partons with small momentum 



fractions |69|| . If shadowing is a small correction, it is a higher-twist effect. Higher-twist need 
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not mean small, however, because when the overlap large, it destroys the incoherence at the 
basis of the partonic formulation, and all twists become equally important. 

The field of small-x dynamics is especially compelling because at very small x, but large 
Q 2 , the density of partons is high, even while the coupling a s (Q 2 ) remains small Such 
a dense, but weakly coupled system promises a new testing ground for field theory, with a 
close relation to the physics of nonabelian plasmas. Even before such an asymptotic condition 
is understood, however, there are many plausiable applications of small-x resummation and 
/^-factorization |nj to current phenomenology, whenever a hard scattering occurs at scale Q 2 , 
with Q 2 <C s. For example, much attention has been given to the production of jets and heavy 
quarks in this regime. Discussions of some of these developments are found elsewhere in this 
volume. 

We have barely scratched the surface of the physics of two-scale problems here. The 
representative examples described above, however, may give a sense of how far it is possible, 
and necessary, to go beyond low-order perturbation theory in QCD. In the next section, we 
shall encounter another extension of the formalism, in which long- and short-distances mix, 
even in cross sections with a single hard scale. 



7 High Orders in Perturbation Theory 

Throughout these lectures, we have used the singularities of perturbation theory as a diagnostic 
for long-distance behavior, and as a guide for organizing the relation of short to long distances. 
In this section, we shall briefly discuss yet another aspect of perturbation theory that gives 
hints of nonperturbative structure, its behavior at high orders [FjJ . 

Recall the relation between the total e + e~ annihilation cross section and the imaginary 
part of the two-current correlation function, eqs. fl75|) and (|76|), 

°e°e- ~ / < I J^J^O) | > , (197) 

with J M an electromagnetic current. As in eq. ( |119| ), we can apply the operator product ex- 
pansion, but now, because there is no "external" momentum in the matrix element, only the 
expectation values of scalar operators can contribute, and only a few operators appear with 
singularities at x 2 = in the operator product expansion, 

! 777 ! 

J"(x) J M (0) ~ -Co(xV)/ + -C q (x 2 ^ 2 )qq(0) + -C F (x 2 fi 2 )F, u F^(0) + ..., (198) 

00 00 00 

with I the identity operator, which does not contribute to DIS. Perturbation theory with all 
masses set to zero contributes at any finite order to Co(x 2 /z 2 ) only. Yet, as we shall now see, 
there is a problem with Cq from high orders, which suggests the presence of the higher terms 
in the operator product expansion, even in the absence of explicit quark masses. 
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Our reasoning begins with the pinch surfaces of Cq, which we have already identified in 



Sec. 3.2, fig. [13], in which a purely soft subdiagram S, consisting entirely of zero-momentum 
lines, is attached to a single hard subdiagram H. Near the pinch surface, all soft momenta 
may be neglected in H, which may therefore be treated effectively as a gauge-invariant local 
vertex. Consider any single gluon internal to subdiagram S, whose momentum we label k. S 
may formally be written as 

s = J d4k j^r t «p ^ > ( 199 ) 

where we have isolated the k propagator (in Feynman gauge) and where T a p is the remainder 
of the subdiagram, as in fig. |28|. Dimensional counting and gauge invariance then require T a p 




Figure 28: Reduced diagram for a soft pinch surface. 

to have the form 

T a p(k, //, a s (//)) = (k 2 G 2 {k 2 )) {k a k p - k 2 g a(3 ) t (k 2 /fi 2 , a s {^ 2 )) , (200) 

where G2{k) is the normalized trace of the full gluon propagator. Neglecting for simplicity any 
renormalization of the vertex H, the function t is renormalization-group invariant, 

d 



ti^/fJ 2 , a s (/U 2 )) = 



(201) 



from which we conclude that we may choose /i 2 = k 2 as we integrate over k, at the cost of 
running the coupling to the scale k 2 , 



(202) 



t(k 2 /fi 2 ,a s (fi 2 )) = t(l,a s (k 2 ) 
= ^2aia s l (k 2 



with the dj numbers. We can already see something funny by looking at the resulting a\ term 
for S. Using the one-loop running coupling, we find 

,a s (k 2 ) 



-3 / d A k k 2 — — 

Jo k z 
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1 , a s (Q 2 ) i , fc2 

1 H — °2 in 7^ 



(/i \ n 1 i 
a s — / dxxln™- 
4irJ Jo x 

= -67rQ 4 a,E ( a ^) T ( n + 1 )' ( 203 ) 

where in the last two lines we have reexpanded the running coupling in terms of a s (Q 2 ) = a s . 
As a result, the nth order in a s has a coefficient that grows like T(n + 1) = n\. This is despite 
the infrared safety of the matrix element. Evidently, this uncontrolled growth in perturbative 
coefficients is a direct reflection of the singularity in the running coupling. All is not lost, 
however, although this behavior will require us to reevaluate how we regard the perturbative 
expansion in QCD. 

A very useful conceptual tool for treating high orders in a perturbative series is the Borel 
transform. Consider a general power series in an expansion variable, in this case a s , 



oo 



n(a s ) = '£c n a a n , (204) 

n=0 

with the c n constants. If H(a s ) is analytic at a s = 0, the c n are coefficients in a Taylor series 
expansion. This need not always be the case, however. We are interested in the case when the 
c n grow, and II may possess no radius of convergence at all about a s = 0. Nevertheless, there 
is a good deal of information in the expansion ( j204| ) . To see why, we define the Borel transform 
of II(a s ) by 

m = e C -T • (205) 



n=0 U - 



It is an expansion in a conjugate variable b, whose expansion coefficients are simply c n /n\. ff (fo) 
is thus much more convergent than n(a s ) is, and has a finite radius of convergence about 6 = 
even when the c n grow as fast as n\. Formally, the inverse transform from if back to II is 

U(a s ) =a~ l H dbe~ b/a °n{b), (206) 



since the integral over b precisely generates n!a™ +1 from the b n term. Factorial growth in the 
expansion coefficients of II now shows up as a singularity in n. If this singularity is on the real 
axis, the inverse transform is an ambiguous integral even if n is known as a function. Even 
if the singularity is off the real axis, its presence indicates contributions to II that cannot be 



described fully by the series ( [204]) , i.e., nonperturbative contributions. 

Returning to our example S^ 1 ' above, we observe that a simple change of variables, 

b' = \n(Q 2 /k 2 ) , (207) 

Air 

in the second line of eq. (|203|) leads to an expression for that is precisely of the inverse 
Borel form, 

r oc p ~87rb'/a s (Q 2 ) 

S^(Q 2 ) = -48ir 2 a; 1 Q i j Q db' & . (208) 
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Here 1/(1 — &2&') plays the role of the Borel transform of S^, and its singularity is a direct 
reflection of the singularity in the perturbative running coupling at k 2 = A 2 . Any such singu- 
larity in the plane of the Borel variable due to the infrared behavior of the running coupling is 
called an infrared renormalon. 



Although ( 208 ) is ill-defined, we have gained something by reexpressing the integral in 
this fashion, if we regard the singularity as an ambiguity in the inverse Borel transform, which 
is well-defined in the full theory, although not in perturbation theory alone. We imagine that 
the Borel transform is well approximated by perturbation theory up to b' = l/6 2 - Although 
the perturbative integral is not well-defined beyond this point, the integrand is already sup- 
pressed at b' = l/b 2 by a factor exp[— S7i/b 2 a s (Q 2 )] = (A 2 /Q 2 ) 2 . Thinking back to the operator 
product expansion, eq. ( 198|) , we recognize that this is the power corresponding to the gluon 
condensate, the vacuum expectation of the operator F 2 . Perturbation theory itself thus sig- 
nals its own incompleteness by generating an infrared renormalon ambiguity at precisely the 
leading nonperturbative power of the operator product expansion (at zero quark mass, in this 
approximation) . 

It is a widely accepted viewpoint that the correct way to treat the perturbative expansion 
is to define perturbation theory by regulating the inverse Borel transform in such a way that it 
introduces a new nonperturbative parameter that may be associated with the vacuum expec- 
tation value (0|F 2 |0). The theory in principle then gives a consistent picture of the function II 
up to corrections of order Q~ 6 relative to the leading power, and up to the next uncalculated 
order in perturbation theory |Tj]] . 

The above discussion has brought us to the threshold of nonperturbative physics, which we 
cannot expect to cross without nonperturbative methods. In closing, we may note that infrared 
renormalons appear not only in the total cross section for e + e~ annihilation, but in many other 
cross sections as well. They are particularly interesting in resummed cross sections, where we 
have seen integrals over running couplings, analogous to those just encountered, appear in the 
organization of large corrections [|72fl . It is natural to ask whether here, as above, perturbation 
theory is signalling a new set of nonperturbative parameters, which probably cannot be reduced 
to the operator product expansion. The full answer to this intruiging question is not, to my 
knowledge, available at present. 
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Appendix A. Color Matrix Identities and Invariants 

I will review in this appendix a few of the group identities useful for elementary perturbative 
calculations in QCD. The "defining" generators are N 2 — 1 N x N traceless hermitian 
matrices, while the generators are defined by the SU (N) structure constants C abC) through 



9309888. 



(T { a A) ) bc = -iC ( 



abc ■ 



(209) 



In representation R, the generators satisfy the trace normalization 



Tr [T^T b {R) } = T R 5 ab , 



(210) 



where for the defining (quark) and adjoint (gluonic) representations 




(211) 



The generators also give rise to invariants Cr, 



n 2 -i 



£ (Ti R) ) 2 = C R I, 



(212) 



a=l 



with I the identity matrix. For the defining and adjoint representations, these are 




(213) 
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For the defining representation, products of the generators in SU(3) are given by 

TPTP = l - [iC abc TP + d abc TP] + l -5 ab I , (214) 
with / the 3x3 identity, and the d abc real and totally symmetric. 



Appendix B. Time Ordered Perturbation Theory, Gener- 
alized Unitarity and the Landau Equations 

Time-ordered perturbation theory (TOPT) allows a simple proof of the generalized unitarity 
discussed in Sec. 3.1. 

Time-ordered ("old-fashioned") perturbation theory is equivalent to the more familiar 
expansion in covariant "Feynman" diagrams. Schematically, for any Green function, 

Green function = E G = E T. (215) 

cov. graphs TO graphs 

Time-ordered (TO) diagrams are topologically identical to covariant diagrams, but their vertices 
are ordered in time. Thus, a covariant diagram with V vertices corresponds to as many as V\ TO 
diagrams. (When a subset of these V\ permulations are identical, they are counted only once 
A TO diagram V consists of the integral of a product of factors, "energy denominators", 
that measure the virtuality of a set of states, 

Here the set of lines between the ath and (a + l)st vertices define the "state" a. E a is the 
"energy of state a" , the total energy that has flowed into the diagram up to the ath vertex. S a 
is the "on-shell" energy of state a, which is the sum of the mass-shell energies of each of the 
lines in a, 



Sa= E <»J= E VlP,l 2 + ^ 2 - (217) 

lines j j in a 

in a 

The factor iV represents "numerator" factors from, for instance, fermion propagators and three- 
gluon vertices, computed with on-shell line momenta. In gauge theories, there are further tech- 
nicalities and modifications associated with extra gauge propagators and self-energy diagrams, 
but we shall not need these subtleties here. The general form of eq. ( |216| ) is, hopefully, familiar 
from TOPT in nonrelativistic quantum mechanics. 



The expression (|216[) is relatively easy to prove directly from covariant perturbation the- 



ory. The example of a scalar self-energy diagram, fig. |29] already illustrates the general pattern. 
In this context, TOPT emerges as the result of carrying out the energy integral(s) of the diagram 
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Figure 29: The two TO diagrams corresponding to the scalar self-energy. 



according to a particular prescription. This consists of (i) reexpressing the energy conservation 
delta function at each vertex as a (time) integral of a phase, (ii) ordering the times of each 
vertex, and (iii) carrying out the energy integrals of each line. For the scalar self-energy, these 
three steps are illustrated by the following: 



I{p,p') 



oo 



dkn I r°° dk' Q 1 



oo 2ir k 2 — k 2 — m 2 + ie J-oo 2tt k' 2 — (p — k) 2 — m 2 + ie 

x (2tt) 2 (5(po - k Q - ko') S(k + k' - p' ) , 
00 dk 1 r 00 dk' Q 1 



2n kl — k 2 — m 2 + ie J-00 2n k' 2 — (p — k) 2 — m 2 + ie 

poo poo 

x / dn e - iTl( - P0 ~ k °-^ / dr 2 e - iT2( - ko+k o-^ 



J —OO J —OO 

1 r dn r dT 2 e- i(p ° Ti - p ° T2) 
j— 00 j— 00 



2uo 2u' 



-00 j —00 

in - T2 )e- l{uj+w ' ){T1 - T2) 



9(t 2 - Tl )e 



-i(w+£j')(T2— n) 



(218) 



where u = ut and u' = uo p -k- The expressions corresponding to the time-ordered scalar 
diagrams shown in fig. result by carrying out the r integrals in (BT 



-i(27r)5(p - p ') [ + l - ] . (219) 

Pq — UJ — to' —UJ — LU ' — Pq 



It is clear that this pattern extends to all orders. 



Within the TOPT formalism, generalized unitarity, represented by fig. [11] and eq. (ff3|), is 
straightforward. To demonstrate fig. O, consider the cuts of a diagram V with A + 1 vertices. 
Applying the energy integration procedure outlined above to the cut diagram, we generate a 
set of cut TOPT diagrams, in each of which a state m is on-shell, with energy denominator 
replaced by 6(E m —S m ), which separates a subdiagram (states 1 to m— 1, denoted r m ) computed 
according to eq. (|216|) and a subdiagram (states m + 1 to A, denoted T^) computed according 
to the complex conjugate rules. The sum over m, for a fixed relative ordering within F is, 
suppressing loop integrals and overall factors, 

A A 1 m— 1 1 

ECr m =E II p — ^—-(27r)5(E m -s m ) n F q , ■ , (220) 

m TO =1 j=m+l ^3 ~ b 3 ~ lt i=l ^>i + ^ 
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where we have suppressed overall factors and integrals. At the same time, the imaginary part 
of (— i times) T, suppressing the same factors and integrals, is 

A 1 A 1 



21m (-iT) = -i 



n 



+ n 



3=1 E i 



le 



(221) 



jJi Ej - Sj + ie 

The expressions ( |220|) and (|221|) are equal, as may easily be verified by repeated use of the 
distribution identity 



1 



2ttS(x) 



(222) 



j + ie x — ie, 

The equality of ( |220|) and ( 221 ) is equivalent to fig. [TT], and holds at the level of the integrands 
of TOPT. The unitarity relation thus holds, as promised, for fixed spatial momentum integrals. 
Only the energies need be integrated. 

New insights into many other theorems of perturbation theory may be found by reconsid- 
ering them in TOPT. An example is the Landau equations fl5"T|). An arbitrary TOPT diagram 
may, following the procedure of eq. ( |218| ) above, be written as an ordered time integral, 

d?L 



r(p) 



dr n ... 



r-.i 



dTo 



d n n 

loops i 



n — 

lines j 



(2tt) 3 r 1L 2u, 

\ ' lines / J 



x exp 



n-1 



(223) 



i (Sai^i) - E a - ie)(r a+l - r a ) - i(E in - E out )r n 

states a=l 

where we have made the ie prescription consistent with Wick rotation explicit, and have ex- 
hibited the loop integrals. The Landau equations emerge as the conditions of stationary phase 
with respect to the loop momentum variables, 

d 



o 



or, 



E E 

i lines j in a 

where v j is the usual relativistic velocity, 



Vj{T a+ l - r a )e\f 



(224) 
(225) 

(226) 



.(a) 



+1 for li flowing in the same sense as the momentum pj in state a; it equals —1 when the 
sense of flow is opposite with pj in state a, and it is zero for pj independent of £i and/or pj not 
in state a. We recognize eq. ( |225| ) as the Landau equations in terms of Feynman parameters 
ccj, by identifying 



E 

a (j in a) 



a 



J > 



(227) 



that is, the ratio of the total time of the states in which a particle propagates to its energy. 
Thus the Coleman-Norton physical process interpretation appears naturally in the context of 
TOPT. 
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Appendix C. The BFKL Kernel and Power Behavior 



In this appendix, I will sketch the derivations of the BFKL kernel in eq. ( [L85|) , and of the power 
u{v) in eq. (|I88|) . 



The vertex in fig. |27] is given by the combination of three-point vertex and eikonal dia- 
grams, 

—gCa'ai 



V a > 



ab\v 



h 2 h' 2 



(k' + k) v g^> + (k + £)^g^ u -{£- k')^^, 



I ]J 2 <lv9m' j. 2 , 



Vv 



W) 



U2 U 2 
rX/rprvjij-t 



U-q U-p 

(k + k) v H p v q v 

s s 



krp 



krp 



Pv 



gC a - 



ab 



U2 y 2 



q ■ £ p ■ £ 

/if I i \ /2£ ' Q 
{k + k) Tu - (- 



V 



Pv 



+ 



,2£-p k' 2 



r 



s q ■ 

where in the second equality we have used q ■ k <^ q ■ £ and p ■ kl -C p ■ £, and in the third 



(228) 



(k + k ) v ~ (k + k )tv H q v Pv 

s s 



(229) 



All of these identities follow from the strong ordering assumption, eq. ( |1 77] ) . 

In the square of the real diagrams, the color factor is easily seen to give — N. The square 
of the momentum factors requires some algebra, and use of 

A£ ■ qi ■ p 



(k - k') 



l\2 



T ■ 



which results in the surprisingly simple relation 

{k' + k) Tv -^-^-)p v + ^-%) qi rl 



s p ■ £ 

The real-gluon contribution to J~{k') is then 



(230) 



(231) 



•F real(^ ) 

w 



-Ng> 



d 4 



(2vr) 4 
a JV fP + dk + 1 



T{k) 



d 2 k 



T 



-F(k) , 



7T 2 Jw+ k+ k' T 2 J (k T - k' T ) 2 kf vv ' ^ 232 ^ 

which, after the redefinition T{k) = !F(k)/k^, and a logarithmic derivative with respect to k' + , 
gives the first ("real gluon") term in the kernel of eq. ( |185| ). 

The virtual gluon contribution comes only from the two-eikonal diagrams of the type 
shown in fig. [27]. In these diagrams, the eikonals play the role of the jets, and the k~ loop 
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integral may be closed on the pole of the lower (p) eikonal to give the logarithmic integral in 
k + . The color factor is again — N. The contribution to T is 



virtual 



(k') 



1 a,N rp + dk+ 



d 2 k T k' T 2 



Hk) 



(233) 



2 7i Jk'+ k+ J (k T - k' T ) 2 k 2 T 

which corresponds to the second term in (|185|) . 

The computation of the power u(v) is reasonably straightforward in dimensional regular- 
ization. Substituting the ansatz (|187|) into the BFKL equation ( |185|) , we get 



Ti 



(k T ' 2 )- iV+l/2 d 2 ~ 2e k T 



- k T ' 2 



1 



d^kn 



{ki)- iv +y 2 {k T - k T ') 2 

1 



k T 2 (kr - k T ') 2 



(234) 



in which e = 2 — n/2 for n dimensions. Both of these integrals have infrared divergences when 
the infrared regularization is removed, but their combination is finite. 



We now use the identity for generalized Feynman parameterization, 

A a Bf> r(a)r(/3) Jo 



(235) 



\yA + (l-y)B]°+f>' 

which holds for complex a and (3. This enables us to perform the transverse integrals by 
standard methods to get 



a s N 
2tt 



nk 



l2\-e 



dyy 



-iv-1/2 



T(-iu + 1/2 + e) 

V(-iv-]/2) Jo -y(\-y)}- % 
dy 



1 f 1 

"9 r ( 1 + e ) / 777 

2 Jo y(l 



ll-e 



The y integrals now give beta functions, in which the infrared poles manifestly cancel, 



(236) 
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r(-ii/ + 1/2 + e) 



B(-e,w + l/2-e) 



7rk 



/2\-t 



r - 



r(-^ + i/2) 

r(-ii/ + l/2 + e) r(iV + l/2-e) 



.-r(l + c)B(-e, -e 



7T 
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(237) 



(The factor (k' 2 )~ e is an artifact of dimensional regularization, which does not contribute to 
the final answer, since in n ^ 4 dimensions the kernel is not dimensionless.) Expanding about 
e = by using 

r(l + 5) = 1 + ^(1) + ... , (238) 
we readily derive the explicit expression for oj(v) in eq. (|188 ). 
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